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Foreword 


Beginning with the issue for July 1959 the Journal of Research of the 
National Bureau of Standards nas been divided into four sections. The great 
diversity of work in progress at the Bureau makes such a step desirable.  See- 
tion B, Mathematics and Mathematical Physics, of which this is the first 
issue, will be devoted to those of the areas indicated by its title in which the 
Bureau is active. These include not only such applied fields as numerical 
analysis and methods of numerical computation, compilation of mathematical 
tables, statistical experiment design, operations research, and the classical 
theories of applied mathematics, but also topics in pure mathematics which 
have a bearing on these fields, like functional analysis, approximation theory, 
combinatorial analysis, theory of matrices, inequalities, theory of games of 
strategy, finite projective geometries, mathematical logic, and the theories of 
stochastic processes and nonparametric estimation, among others. They in- 
clude the theoretical foundations of fluid mechanics, elasticity, celestial me- 
chanics, magneto-hyvdrodyvnamies, statistical mechanics and thermodynamics, 
electromagnetism, optics, and nuclear physics, as well as the underlying math- 
ematical treatment of differential and integral equations and special functions. 
Finally, they include the logical design and programing of digital computers. 

In these and related fields, section B of the Journal will be devoted pri- 
marily to reporting on research at the Bureau. In addition, there will occa- 
sionally be invited papers by authorities outside the Bureau dealing with 
related subjects. We also plan to present occasional review articles, tables of 
useful data, and compilations of information on subjects closely related to the 
Bureau’s program. Thus, section B is planned to include some of the ma- 
terial which heretofore the Bureau has been publishing in its nonperiodical 
series of publications. A complete listing, with selected abstracts, of NBS 
publications and of articles by the Bureau staff in professional journals will 
also be a regular feature. 

C. H. Pace, Editor, 
Section B. Mathematics and Mathematical Physics. 
I 
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Relations Between Summation Methods and Integral 
| Transformations ' 
Werner Greub’ 


(May 12, 1959) 


felations between the Lototsky method of summation and those of Borel, Euler, and 
Knopp are obtained by associating an integral transformation with the series transformation. 


1. Introduction 


In a recent paper, R. P. Agnew [1] ® establishes re- 
lations between the Lototsky method of summation 
for divergent series and the classical summations of 

* Borel, Euler, and Knopp. 

It is the purpose of this paper to show that these 

relations can be obtained in a very natural way if the 

. series transformation is associated with an integral 
transformation. The following scheme deseribes 
the connections between the series and the integral 
transformation: 

original series . > transformed series 


” a 


Laurent expansion Taylor expansion 


+ integral transform 
of the function 


analytic funetion — 


In section 2 it is shown that under certain con- 
ditions a regular transformation can be 
assigned to every integral transformation, and con- 
versely. As an example, the integral transformation 
corresponding to the Lototsky series transformation 
(L-transformation) constructed. This integral 
transformation is used in section 3 for the discussion 
of the relation between the Lototsky and the Borel 
integral methods. The outcome of this discussion is 
that L-summability contains Borel integral-summa- 
bility (Bl-summabilitv) but does not contain a 
modified Borel summability (BI*-summability as 
denoted in [1]). It is shown that the L-method can 
be extended to a summation method which is equiva- 
lent to the BI*-method. 

In section 4 it is proved that a power series is 
L-summable in all interior points of the Borel polygon 
‘of the corresponding analytic function but not 
_ L-summable in all the exterior points. This proof 

is based on the relationship between the L-method 

and the Bl-method. 

In section 5 the Lototsky method is compared 
with the method of Euler-Knopp. 

In the last section the inverse of the L-integral 
transform is constructed and used to obtain’ the 
inverse matrix of the Lototsky series transformation. 


_——— 


series 


“ 


Is 


This work was carried out in part under a National Bureau of Standards 
contract with The American University 
> Present address: University of Ziirich, Ziirich, Switzerland 


Figures in brackets indicate the literature references at the end of this paper. 


2. General Relations Between Series and 
Integral Transforms 


Let a,, (neo=1,2,...) be an infinite matrix whose 
elements are complex and subject to the condition 


lim 4 |a,,,=0 (1) 


We consider a not necessarily convergent infinite 


series 


As 
4 

~ 

~ 


v=1 
with the property 
lim y |u,|< @. (3) 


to « 


This condition assures that the radius of convergence 
of the power series 


is positive. 
Under conditions (1) and (3) all the series 


‘ 


iM , 


U,= 


Gast, 


are convergent, as can readily be seen with the help 
of the root test. Therefore, we can say that under 
those conditions the matrix (a,,) defines the series 
transformation 


(4) 


We will show now that one can associate with this 
series transformation a certain integral transforma- 
tion between the analytic functions 


va = ‘a 
f(z) 2 _ (5) 


and 











‘to obtain the kernel of this transformation we 
define the functions 7’,(2) by 


XN . 
P..(z) ==> Ge 412" (n=1,2.... 
n= 


Because of (1) these series are convergent for every 
zand thus /’,(2) are entire functions of 2. 

It follows from (3) that the function /(2) is regular 
outside of a certain circle, p. Now we integrate 


the product, ?,(2)/(2) along a circle R (R>p), 
and obtain 
oP, (z)f(2z)dz f >) By 94122 
. ar 
2 ri a Qnpty=—2ril’, (6) 


pal 


Under the hypothesis that the function 


is regular in the domain a <1, we derive 


from (6) 


») eS. y ” 
=-7/ l nl . 


n=1 


p K (z,a)f (z)dz 


Hence, defining the function @(a@) for a <1 by 
dla) Pp» U0", 
we obtain the integral transformation * 
O(a) ani $ K(z2,a)f (z dz, (7) 


which maps every function f(z) which is regular 
outside of a certain circle into a function (a) 
which is regular within the unit cirele. 

Conversely, let (7) be an integral transformation 


whose kernel A(2z,@) is regular in the domain 
lel< ©, fal\< land tulfills the condition® A(2,0)—0. 
Substituting for f(z) a function® regular for 2 >1, 
2. U, 
{a= 
r=] « 


we obtain for ¢(@) a funetion which can be expanded 
into a power series within the unit circle, 


_— 
ee a 
n=1 


Q\ a) 


* The path of integration in the 2-plane will always beacircle|z=R (PR 

6 This is not an essential restriction because it can be achieved by replacing 
K(2z,a) by K(z,a)—K(z,0). Then the function ¢(a@) is only changed by a con 
stant 

€ It is assumed that there is no constant term in the expansion; in other words, 


that f(@) =0. 





From (7) we obtain between the coefficients w, and, 
U’, the relations 
’ l * bla) l 2h A Ki(z,a) , 
Ui= x7 P Sori a (351) 4 { auiise f(z)dzde 
I fp p K(2,a) u 
; pe ‘ ‘dzda. (8 


, . , zs . _ 
Hence the integral transformation (7) defines a series 
transformation whose matrix is given by 


(s.:) Pf 


It can be shown that under certain hypotheses about! 
the kernel the series transformation (8) is regular iy 
the following sense: If the series >)u, is convergent,! 


} Ki : outs dzda. 


a 


U 
the series 5007, is Abelian summable (A-summable), 


This means that the limit i 
ee oi 


lim ‘> » a ) 


a~1—0 n 


(a real 
exists and the relation 


lim (SU, a”) >>" 
a >l—0 n=1 r=! 


holds. The hypotheses about the kernel are: 

(i) A(2.e@) is a regular function of = and a@ in thy 
domain 2<©, @ <1 and therefore can be expandes 
in the form 

K(z,a)=>5 a a)z’ 9 
7] 0 


(i) For every fixed n we have 


lim ad,(a@) l la real). 
a >l 0 


(iil) There exists a constant 7 not depending on « 


such that 
¢ Ona) 


—_— 
n=1 


a)—da,(a)| SM. 


To prove the above statement we consider the trans+ 


formation 


Te Ae : 
o(a)=,5 QD K (2,0) f(2)d(2) 
2r. 


and substitute for A(z,@) and f(z) the expansion: 
(9) and (5), respectively. So we obtain the series 


| to-function transformation 


‘ 
— 
t l 


U,,(a). (10 
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Using a well-known regularity theorem we conclude 
from conditions (ii) and (iil) that this transformation 
is regular; this means, the limit lim ¢(@) exists pro- 


a >l 0 
vided that the series S3u, is convergent and the 
T 
relation 
. XX 
lim ¢(a)=>°* u,. (11) 
a »1—0 v=] 


holds. 

The statement which has just been proved can be 
formulated as the following: 

THureoremM 1. Let K(z,a) be a function with the 
prope rties (7), (1), and (717), and let the matrix a,, 


he defined by 


tu (5t;)f PAE det 


Then the series transformation induced by this matrir 
Su, into an A-summable 
t 


series SU’, with the same sum. 


carries every convergent SEriEes 


N 


It is worth noting that the series 5007, need not 


" 
be convergent itself. However, if the matrix a, 
fulfills additional inequalities of the form 
—M 


X 
> [Gal S 
r=1 Mi 


(i - jee ) 


where VJ is a fixed number, it follows that 


—M 


i 


U, 
In this case we can apply the Tauberian theorem and 
conclude that the series 307, is convergent and 


” 


 /, : “ 4 ” 
b Be , lim > (U',a”"). 
7 a>l—0 on 


Thus, under the additional hypothesis (12), we get 


Be i Ss " 
—_— n — 
n c 


for every convergent series >)u,, which means that 


t 
the transformation (4) is regular. 
As an example, let us consider the series transfor- 
mation 


' l ~ 
U,.=wu,, U, . > eo (13) 
@.g¢ 1 


where the p,, are the coefficients of the polynomial 


1)... (z2+n—1) fq= 2... sh. Ce 


This transformation was introduced by Lototsky [3] 
will 


and be called the L-transformation. If the 





series >O/,, is convergent respectively A-summable 
n" 


. . . s 
the original series >3u 


> Su, will be called L-summable 


Ul 
respectively AL-summable, and we write 


>U.=L(S«, 
and | 
/ im ( » » U’ a" ) =AL ( >» u, ) 


For the integral kernel, K(z,a)= >" P,,(2)a", of the 
n=1 
L-transformation, the functions /’,(2) have the fol- 
lowing form: 


? - i-@ 
, i l g (1—a) l 
K (2,0 * (—1)" ( )a"= 

o—1 f= n o—] 

Hence, defining the analytic function (a) 
(@ 1) by 
1 4£(1—a)'-*—1 ' 
o(a)=5 Pp f(z)jdz, (16) 
Zr, z—I 


the AL-summability of the series >)u, means that 


c 
¢(a) approaches a limit if a tends to 1 from the left 
nliong the real axis and we have the relation 


lim dla). 


t a>1—0 


In order to show that a convergent series is also 
AL-summable to the same value we check that the 
kernel A(z, a) satisfies the hypothesis of theorem 1. 
Condition (i) is immediately clear and (ii) follows 
from the fact that @,(@) ean be written as 


1 £K(z,a) ny? 
/ adh, = Y » >. = 
(a1, (a) LL dz=1—é (17) 
ong ag v=1 v! 
where 
v log (l—a). 
To check (iii) we observe that 
yrti 
Ay 4)(a)—a,\a@) e~" ? 
(n+1)! 











ana therefore 


@ x y"rl 
. |dn41(a)—a,(a@)| ” aks p> ’ 
n=1 m=} (? 1)! 
ee" =] (OSa<l). 


Now it follows from theorem | that the transforma- 
tion (13) maps a convergent series into an A-sum- 
mable series with the same sum. 
In this special case it even can be shown that the 
series >35U’, is convergent itself. From (13) we get 
n 


o 1 
Do Envy Dy Pn—1 Pri . 
v=1 Me p=) nh. " 





and thus condition (12) is fulfilled which assures the | 
convergence of the series >07,. This means that 
nN 


the L-method is regular. 


3. Connections Between the Lototsky and 
the Borel Summabilities 


As mentioned before, there is a close relationship 
between the Lototsky and the Borel-integral methods 
of summation. Let Su, be a (not necessarily con- 


t 
vergent) infinite series such that 


n 


li {\Bn!__9 Is 
im «/ : (18) 
" coal n! 
Then the sum 
, | 
, Meas 2 
fi - = { 
F(t)=>)5 xf (19) | 
t=) / 
. . ° : : a : -" . 4 
defines an entire function of ¢. The series S\w, is 


~_— 


called summable with respect to the Borel integral 
method (Bl-summable) if the improper integral 


 ¢e-'F(t)dt 


exists. In this case we write 


BI( Sw.) | ( 


JW 


Fa t dt. 


We are going to show that the function 


ey; 
-| e~ tk (r)d1 
i) 


e 


Bit) (20) 


can be represented as a certain integral transform 
of the function 


~ 
~ 


i 
° 2 
> > ao (21) 


4 


2) that the radius of con., 
; in other words. 


Assume (as we did in see. 
vergence of the series (21) is positive 


that 
—, l 
lim y 4, <o, (29) 
u- r 
; 
Since yn!->o@ it follows from (22) that condition, 


(18) is automatically fulfilled. 


| 


The coefficients u, are represented by the integrals, 


- a 
Uns x flzj2"dz (23 


2ri 


where the integral is taken around certain circle 


} 


R (Rr). From (23) it follows that 
» Mees... l ‘ oa 7 x waew l ‘<a 

ao Mt =5> O f(z) DS | dz=5— — f(zye*dz, 

n=0 ft <7! . a= 0: =-7/. ‘ 

and therefore, using (19) and (20), 5 
' lt ewer } 
Bit) . p @° f(z)dzdr. (24 

~7! e/W e« 





Since the integrand is regular we can interchange the 
two integrations and integrating with respect to 7 we 
obtain 


1 fe@-V1 | ad 
B(t) Dei | : | flz)dz. (25 
Thus the Bl-sum of the series S3u, can be repre- 
| sented as the limit 
B] >.) lj | es l #2): ial 
PE um 5 a H2)dz. (20 


Comparing transformations (16) and (25), we obtain’ 


the following relation between funetions (a) and 
Bit): 
o(a)= B(— log(1-—-a (ai<1), (27) 
which can also be written as 
B(t)}=o(1-—« (28 
From this relation it follows immediately that the 


existence of one of the limits 


lim O(a) or 


lim B(t) 
a >l—v tf» 
and these limits 
This means in terms of series: (nder the 
(22) AL- and the Bl-summability 


implies the existence of the other, 
coincide. 
hypothesis 
equ valent. 
Here the question arises whether this equivalence 
still holds without condition (22). We are going 


to show that the implication BI-+AL is true in 


| general but not conversely. 
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We first prove 
LemMMA 1: Let u, be an infin ite series such that the 


function 


F (thay tet! pe 
preaet v 


is reqular in a certain neighborhood of t=O and define 


the functions B(t) and (a) by 


1 
B(th= | Firj)e-tdr 


and 


B(—log (1—a)). 


(a) 
Then the coe fier nts CU", of the erpansion 
‘UT a* 


4 j n 
i 


o\a) 
are the L-tran sforms of u,. 


Proor: The coefficients (7, are defined by 


Bi > Ua". 


“ 


log (1 a)) 


Differentiating this relation we get 


B’ | log (] a@)) » nl’ a . (29) 
a “SF 
and hence, observing that 
B’(t) Fitye 
and therefore 
B’(—log(1—a)) =(1—a) F(— log (1—a)), 
we conclude that 
, ( )" 
2, Ww at's) , — flog (1 —a)]*uy 4 
n c=0 ( 
Using the residue theorem we conclude that 
p ay A= EY A |log (l—a)]" , > 
U. , (Pp ee da re 
Mt o=0 v a a’ 
which ean be written as 
; ] ~~ l Ss » 
U, G_—| UU,41 > ee (30) 
Me. p=0 nN: » l 
where the coefficients a,, are defined by 
' 2 flee r 
n! og (l—a) 
oe 1)° : p | rs | da: (31) 
v ar 





It has to be shown that the numbers a,, coincide 
with the numbers p,, defined by (14). First of all 
we have a,,=—0 for v>n since the functions 


flog (1—a)]" 


a” tl 


are regular at a=0 provided that r2n+1. We now 


define the polynomials a,(2) by 


n 


@,(2) =D) Bae2" 
t 
and show that 
a,(—k)=0 (k=0,1,...."—1). 


In fact we have for every integer 


a,(—k) > 2 fp llog G—wf da 
Nn a 
n! =i v! : att! 


a A. |log (1 a)|" . (l—a)*¥—1 

1 WY —— da=Q : da. 
—_— ‘ a n+l 

gay 2. a Ki a 


If now ks n—1 there is no residue at a=0 and we get 


(k (32) 


QO,l,....m—1). 


From here it follows, since the polynomial @,(2) has 
the degree n, 
(33) 


1)... (2+n-—1) 


where \ is a constant. Substituting z=1 in (33) we 
get 


a,(1)=An! 


and on the other hand we have 


n e >. (—1)" £ [log (1—a)]* 
X Xx ' XN os = . 

a, (1) => ay, => On =2! > fp = da 

e=) r=1 v=1 ve. a a 

y ] da 
mn ( — )- i =n! 
l—a a 

From the last two relations it follows that \=1 and 


we obtain 


1)... (2-+n—1) =p, (2) 


and comparing the coefficients 


On, 


Poo» 





which proves lemma 1. 








THEOREM 2.7 A Bl-summatkl series is also AL- | 
summable and the two sums coincide: 


BI(Su,)=AL( Su.) 


t ys 
Proor: By hypothesis the series 


tl pe 


Sy U, 
— yl 


F(t) 


r=0 


defines an entire function. Hence the funetion 


Bit) | Firje-"adt 


i] 
is entire too, and therefore the function 


a)) 34) 


ola) 


B(- log (1 


According to lemma | we have 


and the sequence uw, by 
Uy, yg” 1 (()) (Vi hb - (36) 


Then we have 


wy Urs ie | 

ra J! ) 

r=0 e. és ( 
and therefore this series is convergent only for 
t<log 2. Hence the Bl-method cannot be applied. 


To obtain the L-transform of sequence (36) we 

A hl a bad > ™ . 

apply lemma 1.) The function (ft) in our case is 
given by 





is regular for ja <1. 
the expansion 


2 


Ss. Ua" (335) 


(a) 


n=l | 


where the l’, are the L-transforms of the u,. Hence | 
the series | 
S30 | 

n@ | 

n=1 | 

is convergent for a@)<1. Furthermore, we obtain | 


from (34) and (35) 


B(—log (1 


> / a" 


n=l 


—a)), 
and hence, if a tends to 1 on the real axis. 


lim Bit). 


to+o 


lim (> l "a" ) 


a>l—" 


This means that the series SOw, is AL-summable 


and that the relation 


AL(> uy ) -BI( 


t 


holds. 


Conversely, one cannot conclude that an AL- 
summable (even an L-summable) is also 
Bl-summable, as the following example shows: 


series Is 


Let the function g(t) be defined by 


l 
g(t=5 , 


- ( 


? Different proofs of theorems 2 to 4 are given in a mimeographed report by | 
R. P. Agnew (U.S.A.F. Contract 18 (600-685) ). | 


; rf 
e ‘dr , 
Bit) , =log(2—e-'), 
02 f . 
and from here we get 
o(a) = B(—log(1—a)) =log(1+a), 
and thus 
l 
if | lies 
i 


The series aul 
2 
is L-summable to the value log 2. 
As the above example shows, the AL-summability 
of a series S3u, does not imply that the series 


a convergent, and hence series (36) 


has an infinite radius of convergence. It can be 
shown however that this radius can not be less than 
log 2. To this purpose we prove 

LEMMA 2: Let ou he an infin ite series and ae "ats 


’ N 


L-transform. Assume that the series 


d(a)= DSS ,a 
is CONV rgent ina certain circle a p, (p >) Then 
the series 
’ ul 
y Ss ' 
I t — yl f 
is convergent in a certain circle t <p’, (p’ >0) and the 
relation 
¥ Urs os ; 
| (x : rye dr=o(1—e!) 37) 
0 vr 
holds. 
Proor: By hypothesis the function 
Bt) =o(1—e' (38) 


is regular in a certain circle ¢ < p’ and therefore the 


same is true for the function 
(39) 


F(t)= B’ (the' 


5 


whe 


Fro 


Exp 


we ¢ 


and 


asa 


On 


Fro 


and 


U) .= 


ny 
tral 
han 
Ur. 


vr 
The 


and 
ere 
R 
(3S 
I 
whi 
1 


ser i 


T° 


where B’ indicates the derivative of B. 


From (38) and (39) it follows that 


| 





Fit)=¢@'(1 e) fj (40) 
Expanding the function F(t) into a power series, 
’ X 
raj=-> ZF, 
¥ 
we obtain from (40) 
@d’ (1—e-*) oe. 
0 
and if we introduce 
a=1—e-! 
as a new variable, 
¢’ (a) =>35(—1)*e,[log(1 —a)}’. (41) 
t 
On the other hand we have 
od’ (a) > nl ,a"'. (42) 
n 
From (41) and (42) it follows that 
>> wl’ ,a"-'=D5 (—1)"e, [log (1—a)]’, 
a ti 
and from here we obtain 
, 1 | 4. |log (1—a)]* 
Runs fe a’ 
_| Ss , 
ee oe Pn pele 
Mt. r 
This equation states that the (, are the L- 
transforms of the sequence vile, On the other 
hand, they are by definition the L-transforms of 
u,, and hence we conclude 
l 
¢, u, 
vr ' 
Therefore the expansion of F(t) becomes 
. my Ue} ‘ 
F(¢)=—2° = t* (43) 
t 0 . 
and this series must be convergent within the 
circle ft). p’. 
Relation (37) now follows from (43), (39), and | 
(38). 


From lemma 2 we can derive the following theorem, 
which is proved by Agnew [2] in a different way. 
Tureorem. For the terms of an AL-summabl 
series >3u, the inequality 
’ 


lim “ (44) 
n> @ \ ni log » 


TOSYS6 5 2 7 


holds. 


Proor: Since the series >}u, is AL-summable the 
function ’ 


Qa) 2 


" 


is regular for |a <1. Hence the function B(?) 
defined by 
B(t)=o(1—e™) 


must be regular in the domain LD), which is deter- 
mined by the inequality 


(45) 
and the same holds for the function 

F(t)=B’ (t)e*. 
Writing t=r-+-71y, we get 


|1—e-'|?=1+4 2e-? (} e~7—cos y) 


and hence inequality (45) is equivalent to 


Se77<cos y. 
This means that the point t=r+7y lies in the 
“interior” of the curve 
r=—log (2 cos y). (46) 


y 
a 
— 


—> xX 





_——’ 
log2 





Pica 


Relations between summation methods and integral 
transformations. 


Figure 1. 


Thus the function B(t) can be expanded into a 
power series within the largest circle f= p with the 
center ¢=0 that is contained in the domain D. 
This is the circle 

\t}=log 2. 
because this circle touches curve (46) in the point 
t log 2, it is sufficient to show that the radius of 
curvature of the curve (46) is not less than log 2. 
The curvature of this curve is given by 


l l 


p cosy 











and hence we have 
p>1>log 2. 
Thus the expansion of F(¢) into a power series 
must be convergent for 


|t|< log 2. 


On the other hand, using lemma 2 this expansion is 
given by 


a ‘ 1 : 
F(t) > & r 
, 0: 
From here it follows that 
—_—— © uy, 
lim | + . , 
we V nt “log2 


which proves theorem 3. 
It follows from series (36) that the bound 1/log 2 can 
not be improved. 








The BI*-method. The application of the BI- 

method provides that the series | 

“ U, +1 

>> 7 © 

; & 
. cine % . —— . 
converges in the whole f-plane. Frequently the | 
case occurs that the radius of convergence of this | 


series is positive and finite such that this series de- 
fines a regular function F(t) within a certain circle 
tip. Ifthis function can be analytically continued 
along the positive f-axis one can form the improper 
integral 


| F(tye-'dt. 
J 0 


In case this integral is convergent, the series is said 
to beBL*-summable [2, sec. 11] and its value is called 


the BI*-sum of >}u,, 


BI*(S>u,)= [Fin dt. 


Obviously a Bl-summable series is also BI*-sum- 
mable to the same value. 





It can be shown that AL-summability implies 
Bl*-summability, but we do not prove this here 
since BL*-summability will turn out to be equivalent 
to a generalization of AL-summability and the above 
statement will follow immediately from there. 

To demonstrate the BI*-method by an example, 
let u, be defined by 


' 
ig 


Un+1 = an 


° . } 
where a@ (a0) is an arbitrary complex number, 
In this case we obtain ; 


and thus, 


BI* (>) | 


aft ; 


This integral converges for every number a which ! 
is not real and negative and therefore series (47) jis; 
BI*-summable for all those values of a. 

The next example shows that a BI*-summable 
series is not necessarily AL-summable (and hence 
not L-summable). Let the function ¢(¢) be defined by 


olf) ——y 
lad ( )~ 


(48 


where a is a real number such that a@>1, and lets 


the sequence wu, be defined by 


' 
u,=o'"~"(0) (se=—12. (49) 

Then the series 

>" . f* 

homed ' 

, ’ 1)! 
is convergent in the circle 

t\< log a 


and represents the function $(f).~ This function is 
defined along the whole positive f-axis and the integral 


is convergent. Thus the series (49) is BI*-sum- 


mable. 
In order to obtain the L-transform of (49) accord- 
ing to lemma 2 we have to form the function 


' l od : 


a—l a—e 


Bit) d(rie~tdr 


and then substitute 


This gives 


QO\a) — ° 
a—1 a-—l-ea 


The function ¢(a@) so obtained has a pole at a=1—a 
and hence the series>)/7,a" is only convergent in 
the cirele ° 


whic! 
fore, 
Th 
the 
meth 
to a 
exter 


Le 


series 


is col 
the 

along 
limit 


exists 


write 


AL*- 


ity. 
TH 

to the 
Pr 


mabl. 


is reg 


tions 


and 


Then 
there 


Acco) 
L-tra 
Fit) « 
t-aXxis 
Bit). 
tion « 
and y 


Ata 


- 


il 


which is smaller than the unit cirele if a<2. There- 
fore, the AL-method can not be applied to series (49) 

The AL*-summability. As the last example shows 
the Bl*-method is more powerful than the AL- 
method. However the AL-method can be extended 
to a method AL* in a similar way as BI has been 
extended to BI*. 


Let Su, be an infinite series such that the power 
t 
series 
g(a) =>)’, a" 


n 


p and assume that 
continued 
the 


is convergent within a cirele ‘a; 
the function ¢(@) can be analytically 
along the interval OSu<1 (a=u-+ir). If 
limit 

lim o\a) (a real) 


a>!l—v 


exists we shall eall the series S‘u, AL*-summable and 


eed 
7 
write 
AL*( ot, ) - lim (a). 
t a >l—v 


AL*-summability obviously implies AL-summabil- 
ity. 


THeoreM 4.) The AL*-summability is equivalent 
to the B1*-summability. 


Proor: Assume first the series SSu, is Bl*-sum- 


“ 


mable. Then the function 
. « Mes 
I (ft) iy . ! t 
- ! 


is regular within a cirele |¢}<p. Let now the fune- 
tions Bit) and o(a) be given by 


alt | 
Bit) Fi rje-tdr 
Jt 
and 
¢(a) = B(— log (1—a)). (50) 
Then @(a@) is regular within a circle ja <p’ and 
therefore it can be expanded into a power series 
d(a)= >)’, a". 
n 
According to lemma 1 the coefficients Ul’, are the 


L-transforms of the u,. By hypothesis, the function 
Fit) can be analytically continued along the positive 
f-axis and therefore the same is true for the function 
Bit). Hence eq (50) defines an analytic continua- 
tion of the function ¢(@) along the real unit interval 
and we have 


lim (a) 


a ?>l—0 


lim B(f). 
lta 





This means that the series >)u, is AL*-summable 
- 


and the relation 
AL*( Su, =BI*( Du, ) 
holds. ' 


On the other hand, if the series S}u, is AL*-sum- 
. 
mable, the function 


d(a)=>) Ua" 
n 


is regular within a cirele |ai<p. Hence the series 


, Ue 
F(t)=)> 4 
vp! 


t 


according to lemma 2 is convergent within a circle 
t<p’, and the relation 

F(t)=9' (1—e~*)e! (51) 
holds. Since the function ¢(@) by hypothesis can be 
analytically continued along the real unit interval 
the same holds for F(t) along the positive t-axis and 
relation (51) holds for every positive ¢t. From there 
we obtain 


Fit)e~‘dt= lim ¢(a) 


J a1—0 


which shows that the series >}u, is BI*-summable. 


D 
Altogether we have the following relations between 
the Lototsky and the Borel methods, where all the 
simple arrows are not invertable: 


L AL 


>A L* 
a» IT 
BI BI* 


Series (36) shows that L does not imply BI, but it 
is not known to the author whether the implication 


BI-+L holds. 
3. Application to Analytic Continuation 


Let g(z) be an analytic function regular in a do- 
main J) which contains the origin and let 


g(z)= 2D) 
v=1 


as" (52) 


be the expansion of g(z) around the point z=0. It 
is a known result that the series (52) is Bl-summable 
to the value g(z) at every interior point of the Borel 
polygon mw of g(z) and not Bl-summable at every 
exterior point. Making use of the relationship 
between the BI- and the L-summability the same 
can be shown for the L-summability of series (52). 
In order to do so we first observe that 
lim ‘y\u,2%-'|=|zollim 4 \u,|< @, 


ui " 











and hence the BIl-summability is equivalent to AL- 
summability (sec. 3). From this remark and the 
implication L->AL it follows that series (52) can 
not be L-summabie at an exterior point of 7. 
Let now 2 be an interior point of z Then the 
series >) u,2)7' is Bl-summable and hence AlL-sum- 
t 


mable. Thus the relation 
lim (20 0,(0)a" ) (Zo) 533) 
al-—0\ Hn 


holds where @7,(2)) denotes the L-transform of the 
sequence U2)” °. 
It remains to show that the series 5007,,(29) is con- 
a 
vergent; then the relation 
i 
230. 


(Zo) =o) 


a 


will follow from (53) and the Abelian limit theorem. 

We are going to prove the absolute convergence 
of the above series. For that we start out from the 
relations 


” 1 
a a ; Ss 
U",, (20) -> 1 Pa-i » Un+i Zo- (54) 
Mey =] 
Since 2 is an interior point of the polygon 7 the 
closed circle 
: ~ < <0 
zi. 2 
is contained in J), and therefore a number p>1 can 


be chosen so close to 1 that the circle 


still is contained in D. (2, Chap. VIIL.| 
Hence the coefficients u, can be represented by 


1 fg(z) 7 
ong dz. (5: 


From (54) and (55) we obtain 


v_! ‘ 92) a, 1 gg) 29 
c 


~t 
~ 


u, 


and irom there, if \/ denotes the maximum of 
(z) . 
gi*) on the cirele e, 
2 
; M Zo 
< 

If we write for abbreviation 

w —y 


| 


10 


this relation becomes 


M 
[ Nn < 9 |@||@ l Ww i 2 
NM. 
Wi Ww | s (Ww Ni) Hv 
M ’ ae * Lae 
nen wn l)(wt+n) 


(56) > 


If n tends to o the first factor approaches I'(w 
and therefore stavs bounded. The same is true 
for the second term and = therefore is sufficient 
to show that the series 


Is 


it 


yy. . y. ~- 
" (4 


7 


circle c 


is uniformly convergent on ec. The is 
represented by the equation 
_ ' 
-O | | a | < <9 
“i = 1 -- pe QOs¢dS2r t 
From here it follows that 
: 9 
1+ pe’? 
and for the real part 
2(1-+ p cos 6) 2 
R(w) “ <s— | 
1+2p cos 6+p*~ 1+ p 
Therefore we have 
») 
: 2p . , 
Riw 2) 3) —1—6, 6 >O) 
and hence the inequalities 
it ‘1sn 


hold, which umply the uniform convergence of series . 
(57) on ce. Now the absolute convergence of the 
series > » follows frem (56). 


/ 


5. Connections With the Euler-Knopp 
Method 


The Euler-Knopp transform of a sequence 4, 
(EK-transform) is defined to be the sequence 
[ ; S % ” c n—t “© 
n D(. a (l—« U, (58) 
t 


where « is a complex parameter. 
This transformation can be derived from a linear 
transformation of the w-plane which has the fixed 


poin 
gene 


whe 


be ¢ 


This 


is re 
und 


Sine 


func 


The 


If 


eXpl 


This 
thec 


whe 


Fre »] 


point w= 1 and maps w=0 into w=. The most 
general transformation of this kind is given by 
1 
a . 
+a (59) 
Ww 
where a is a complex number. Now let 
./.. , U, : 
f(z)=35 - (60) 
’ S 
be an analytic function which is regular for | z # 
This implies that the function 
, .{l—a 
T*(@) f( +a) (61) 
Ww 
is regular in the domain that is mapped onto (2 >1 


under the transformation (59), that is the circle 


Since this domain contains the point w=0 the 
function f*(@) can be expanded into a power series 


f*(w) Se git". 


" 


The coefficients Ul", are given by 


, | < f *(q) 
{ . D dw. 
’ ») n+1 
<7 Ww 


If we substitute here for the function f*(@) the 
expression (61) and for f(z) the expansion (60) we get 


Ww 


: l SO c w \° dw 
U.=< x P( ) —— Up. 
k l—a+aw 


This integral can be evaluated by the residue 
theorem. First we observe that 
, a ; 
(l—a-+-aw) (l—a) (1+ ) 
l—a 
x"; 1 (K lw : 
where 
l 
A= (G2) 
l—a 
From there we obtain 
" —* . c 1 
l =5—-(D > «*( ) 1 )4w?* + do) u 
-Ti “ ron be 


“ 
) a(n iio 
- 
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and making use of the identity 


eee 


we finally get 


I > (" 1? 


) cil—2z)*"" a. 
v—l 


This equation shows that the coefficients U", are the 
KK-transforms of the coefficients u,. 

From the relationship between the EK-transforma- 
tion and the linear mapping (59) it can easily be 
derived that the EK-transformation cannot be 
regular unless « is real and O0<«S1. Assume that 
transformation (58) is regular and let S3u, be a con- 

t 


vergent series such that the function 


f(z)=>5 — 63) 
w ~? 


2,41). 
» by hypothesis is also convergent 


has a singularity at a given point 2, (2 a 
= _ 1 I 
Then the series Sv 


and hence the function 
. jl--a 
fw =S( +a ) (64) 
w 


is regular for jw@<1. Therefore the funetion f(z) 
must be regular in the image domain of the cirele 


w<1 under the mapping (59). This image 1s 
given by 
D,: |2—a\>!1—al. 


Therefore the point 2; cannot be contained in D, 
and since 2, was arbitrary on the unit cirele this has 
to hold for the whole ecirele |(z 1. Hence the unit 
circle has to be contained in the cirele 


and this is possible only if a is real and a<0, which 
means that « is real and 0< «S11, 

It is known that an EK-summable series is also 
AL-summable provided that « is real and positive 
(2, see. 8.) We shall give here a proof of this fact 
using the integral representation of the L-trans- 
formation. 


Let Su, be an EK-summable series and define the 


function f(z) by 











Then the function 
* 4 
g(w) =f ( ) t((1—a)w+a) 
WwW 
has the expansion 


g(w) = >) Uw". 
n 


The function f(z) can be expressed by g(@) as 


per=o(§=8) 


From this equation it follows that the L-integral 
transform of f(z) can be written as 


2? £(i—«)'"*—1, 
o(a)=5— fp = f(z)dz 
2ri c—1 : 


l . (l—a)'-*—1 2—a 
oF” se sai fe) 


— 
+ 


Introducing the new variable 


we obtain 
] . (l1—a) O-Pa-o_ | 
o(a)=5 —s fp - gi Ode 
2m f— 
and hence ¢(a@) can be considered as the transform of 


g(w) with respect to the kernel 


(l—a)"-P-e l 


K( [a)= ™ 
c—1 


"t is easily verified that this kernel fulfills the condi- | 


tions of theorem | provided that ais real and 0 <a< 1. 
Therefore it follows that under this condition 


lim ¢(a)=3> U,=EK (Sou, ) 
a>1l—0 u t 


provided that the series S30’, is convergent. But 
n 

that means that the series >\w, is AL-summable and 

the relation 


AL (Su, )=EK (Su, ) 


. 8 


holds. 

In case «x is not both real and positive it is shown 
by Agnew ([2], sec. 8) that EK-summability does 
not imply AL-summability. , 


6. Inversion Formula of the Lototsky 
Integral Transformation 


The question arises whether the integral transfor. 
mation 





l ~ (l—a@)'-?7—-1, ' 
d(a) =; fp —s fle)dz (65 


given regular function @(@) (a@<1) there exists a 
function f(z) regular outside of a sufficiently large 
circle such that relation (65) holds. Provided such 
a function exists, it can be represented as an integral 
transform of ¢(a@) at least in a half-plane. 


2 >r such that f(@ )=—0, and let d(a) be its transform 
with respect to (65). Then in the half-plane 2 >r th 
function f(2) can be represented as 


“1 
f(z) (z | (a) (1—a)*~*da. 


as 


l 4 (] a) — 
(a) 5-4 a P(z)dz. (06 


Let 2=s9+- iyo be a fixed point in the half-plane 
a circle (2 =R, where 


. R. Jt. 


respect to a between 0 and tf (O0<t< 1), we get 


] 


; “ I **§ £ (l—a)*e-? ; 
(a) (1—a)**" “da —s | | t2jdzda 


=-T7! Jo. 





v7! 


| Since the integrand is regular, the two integrations 
/may be interchanged and we obtain integrating 
with respect to a 


% _ l £(1—t)*e-*—1 , 
| dla)(l a)70-* da — p f(z)dz. 
| Jo omuy (2—1)(Zo—Z) ° 
If ¢ tends to 1, the expression 
(l—f 
tends to zero uniformly on > =R beeause of the 
inequality 
Re (2.—2) =m—R>O0. 
| Hence we obtain 
ae | . . 
l f(z)dz = 
o(a)(l1—a)*0-* da q ° ; (O74) 
J 0 2T i, {2 | )( 2, S 
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) 


ean be inverted. One cannot expect that to every | 


THeorem 5. Let f(z) be a funetion regular for. 


Proor: Since f(@ )=0, relation (65) can be written’ 


s>r. We choose for the path of integration in (65)° 


If (65) is multiplied by (1 —a@)*~* and integrated with » 


The 1 
by tl 


f(2) 


we gt 


Fron 


and 


whic 

I 
mati 
part 


rm: 
rhis 


the 
is re 


we 


The 


The 


and 


ry 


a 
ze 
ch 


‘al 


' 


or 


mh 


he 


NT 


cr 


The right-hand side of this equation can be evaluated 
by the Cauchy integral formula since the function 


f(z) is regular for (2 >r. Observing that f(@)=0, 
we gel 
& f(z)dz f (Zo) 
f : =! : (68) 
2muy (z—1)(Zo—2z) 2-1 
From (67) and (68) it follows that 
f (29) - oth 
Z = | o(a)(1—a)" “da 
<0 l J 0 
and we obtain the formula 
*1 
f(z)=(2z-1) | (a) (1—a)*~*da (r>r) (69) 
J 0 


which is the inversion of (65). 
Formula (69) can be used to obtain the inverse 


matrix of the L-transformation (13). For the 
particular functions 

o(a) =a (k=0,1.. .) (70) 
we get from (69) 

*1 
f(z)=(z2—1) | a‘ (l—a)*~*da 

k! k! 
-2(2+1) (c+k 1) plz) 


This function is in fact regular outside of the circle 
k—1 and thus for the special functions (70) 
the integral equation (65) does have a solution that 
is regula r outside of a certian circle. 
To obtain the inverse of the L-transformation 


(71) 


we denote by g,, the inverse of the matrix p,,,. 
Then it follows from (71) that 


", p> Idn-1 5 C7 (72) 
The funetions 
P\a) ee ia, 
" 

and 

. u 

(2) ie “We 

. oad t ? 

7 
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correspond to each other by the transformations 
(65) and (68). 
In particular, the function (70) 


coefficients 


the 


possesses 


l nN Onks 
and hence we obtain from (72) 


a A 
Uy» => T'Un 1j 15 x kqy- 1k—1- 
J 


Thus the function f(z) has the expansion 


-— ¢ n Lend 
f(z)=ht > et. (73) 
On the other hand this function is given by 
k! . 
f(zaj= (74) 
Py(2) 
Comparing the right-hand sides of (73) and (74) 
we conclude that 
> 2 pe-1__ l 
" a Px(2) 
and from there 
l f ge" 
( sa——; Ui) dz 
ibaa 2ri J (2) 
or if n—1 is replaced by n and k—1 by &, 
l fp 2" ie 
WT» —y dz. (79) 
on 2ri J pai (2) 


This formula corresponds to the representation (31) 
of the matrix p,,. If the integral in (75) is evaluated 
by the residue theorem we obtain the formula 


k pn 1 
q, > (—1)"-° 
Dak ao (yvn—1)'(k—p)! 
l 2. k—1 
2 (—1)* “( yom. 
(k—1)! 2 v—1 
7. References 
fi] R. P. Agnew, The Lototsky method for evaluation of 
series, Mich. Math. J. 4 (1957). 
(2) G. H. Hardy, Divergent series (Oxford, England, 1949). 
(3) A. V. Lototsky, On a linear transformation of sequences 
and series; Tvanov. Gos. Ped. Inst. Ue. Fiz.-Mat. 


Nauki 4, 61 (1953) (in Russian). 
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On a Modification of Watson's Lemma 
F. Oberhettinger* 


May 27, 1959) 


The method of steepest descents is extended to the case when a saddle point and a 


pole of arbitrary order are involved. 
is demonstrated. 


l. Introduction 


The solution of numerous problems in electro- 
magnetic and acoustic wave propagation theory is 
frequently given by an integral representation of 
the form 


g(z) | wire dr, (1) 
Jc 


Here g(z) denotes the desired property (for instance 
the electromagnetic or acoustic potential). The in- 
tegration in (1) is taken along a curve C in the 
complex ¢ plane along which w(v) and u(v) are 
defined. These function contain generally certain 
parameters which depend on the physical configura- 
tion under investigation. The parameter z is usually 
related to the distance between some point of refer- 
ence and the point of observation. A reduction of 
the integral (1) to an asymptotic series has to be 
carried out When approximate expressions for large 
>are desired. The frequently used method of steep- 
est descents [3, p. 86]' becomes invalid when one of 
the above mentioned parameters assumes a value 
for which a pole of w(v) coincides with a saddle 
point of u(r). Methods to cope with such 
have been worked out and discussed [1], [8], [9]. 
The former two give expansions which proceed in 
terms of hypergeometric functions besides a “leading 
term”, while the latter consists of an error integral 
as leading term and an asymptotic series with m- 
verse powers of 2. [t is shown here that the method 
given in [9], which involves a pole of the first order 
and a branch point of order 2 can be extended to 
the case involving a pole and «a branch point of 
arbitrary order. The expansion for this case con- 
sists of a finite number (equal to the order of the 
pole) of terms involving Whittaker’s functions and 
an asymptotic series with mverse powers of 2. 


causes 


2. A Modification of Watson's Lemma 


It is assumed that the reduction of the contour 
imtegral (1) (for instance by the method of steepest 
descents) to one or several integrals of the Laplace 
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An application to a problem in diffraction theory 


transform type 


f(z) -| Fitje-"dt (2) 
J 0 
has already been carried out and only this kind of 


integrals shall be considered. A theorem commonly 
referred to as Watson’s lemma [2, p. 218] admits the 


asymptotic representation of f(z) for large z when 
the behavior of F(t) for small t is known. In its 


general form given by Doetsch [3, p. 48] and [4] one has 
the 

THeorem: Let F(t) be analytic in a sector a< 
arg t<8 of the compler t plane with the possible ex- 


ception of the origin, and let the behavior of F(t) in 
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this sector be 


(a) F(t)|\< Ae*'"', for \t) >R, A,c>0, 

(b) F(t) = S3a,te, —1< Re ap Re a <: - + 
n=0 

when t->0 in ax arg (<p. 


The asymptotic expansion for f(z) = | F(t)e-*dt 
J0 


for large 2 is then obtained by inserting the series (b) 
for F(t) into (2) and integrating term by term. 


f(z) F(t)e—“dt = >)a,T (1+a,)27%, 7! (3) 
eV n=0 
—- = 
when 2-> ©in—5—B< arg 2< 5—a. 


The funetion F(t) in (2) is now assumed to be of 
the form 

F(t) teq(t), Rea>—1, (4) 

where the function g(f) is analytic in the vicinity of 


‘—0. In this case 


f(z | teq(the “dt (5) 


1 fd" (¢ 
a,= sro ( ) ° 
n! dt : g t=O 


with 


g(t)= >) a,t", 


” u 


(6) 








to 








This series is convergent for ft < ft), where fy is that 
singular point of g(t) which is closest to the origin. 
The other singular points are denoted by 4, fy, 


such that their moduli form a nondecreasing se- | 


quence. 
is on the positive real f-axis. The asymptotic ex- 
pansion of (5) then becomes by (3) and (6) 


fi2)= | taq(tye 
0 


for 2 The range of validity for arg z depends on 
the behavior of g(t) according to (a) and (b) of the 
theorem quoted before. The case in which the 
singular point f) is a pole of order m, say, is considered 
now. The Laurent expansion of g(f) in the viemity 
of ty) is then 


—"dt = > a,0(at+n+1)2-2-""'! = (7) 


n=0 


>co,. 


q(t)- b_,.(t—to)-"+ ... +6_,(t—to)-'+9*(t) (8) 
where 
g*(t)=9(t)—b_,(t—to)-"— .. . —b_,(t—ty)=' (9) 


and consequently can be 


d” " 
(7 ’ t ft, - 
E J | 0 , 


(10) 


is analytic for f< ft, 
represented by the series 


2 l 


g*(t)=>) ek", 


n=0 


It has to be remembered that in physical applications 
the function g(t) in (4) contains generally certain 
parameters which depend on the configuration of the 


We assume that none of the singular points | 


and (b) of the theorem quoted before. 


problem, and this may become such that the pole 


t, of g(t) coincides with the origin f=0 (see section 3). 
In this case all the coefficients a, tend to infinity by 
(6) and (8) and the asymptotic series in (7) becomes 
useless although the integral in (7) may exist. It is 
therefore necessary to obtain an expansion which for 


tj #0 represents the integral (7) asymptotically for | 


large z and which for 4;=0 (if the integral exists) 
gives its correct value. Such an expansion can be 
obtained as follows: One inserts the representation 
(8) for g(t) into (7) and gets 
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f(z) teq(the-“dt 


7h. | (t—t) ~*t@e- "dt + | tag*(the-"dt. 
k=1 J 0 0 
(1] 
But 
t*e-"= (t—t))-*dt=T' (1 +a) (—tp) 22734228 22 
ly , 
( ba || 1 ee eae ee faz) (12 


by virtue of the integral representation for Whittak- 
er’s function [5, p. 274]. The integral 


J, 


can be represented asymptotically for large z by 
means of the power series expansion (10) (which is 
valid for |¢/< ¢,!) and Watson’s lemma. 
then finally 


teg*(t)e dt 


f(z) | tey(t)e- “dt 
e o 
, 7 . la Lay 
=> 5 ¢,T(at+n+1)z-¢ '4+ T(1+a)(—t)2%27 247! 
v=0 
m l 
‘ 2 
é . h «( ) WW 1 Sees | >fy) 
I 1 t, 
(133) 


The range of validity for arg z depends on the 
behavior of g*(t) according to the requirements (a) 
For certain 
combinations of @ and m (13) can reduce to a simpler 
expression, as for instance 


‘2tDs_ ,[(22 


(parabolic cevlinder function) — (14) 
Woy y(z)—atetel Enfe (2h) shel [mht (15) 
Wael Le} Kil ty! | te" (16) 
of sae : 
W,..(2) = r)*K,(5 2 ): modified Bessel function). 


(17) 
The case a 5, m—=1 (branch point of order 2 and 
pole fy of the order 1) which leads to the special case 
(15) has been treated before [9]. For this case one 
obtains from (13) and (15) 


One obtains ' 


de 


or 


Bi 


ap 
sh 
di 


go 


Wl 


q| 


la 
T 


of 
of 
te 
m 
pi 


ak- 


by 
1 is 
Ins 


he 
a) 
in 
ler 


4) 


Hy) 


i “— , l n 
f(z)= t iy(tdt~ de (n+5) 2 2 
+-b_yw(—ty)~ 4e7*0 Erfe [(—2ty) 2]. (18) 
Te coefficients ¢, according to (9) and (10) are 
defined by 
q(t) b 1(t—fo) > et © A (19) 
n=0 
or ] 
] ( vw 
. (t)—b_,(#¢—t,)~' . (20 
“a ae ly te }, 0 


But these coefficients remain finite when t)—>0. 


3. Application to the Wedge Diffraction 
Problem 


The investigations outlined before will now be 
applied to a previously discussed problem. It was 
shown [6, 7] that the solution of the problem of the 
diffraction of plane waves by a wedge or corner is 
governed by the integral 


I (a,b) - t-hy(t)e "dt (21) 


i) 


with 
(2a)~'e~ 7? sin (mp/a)(t4-2)~? 


(1?+4-2t)4)]|—cos(mp/a)}~'. (22) 


qit) 


- cosh [t/a log (14-44 


Here g(t) is analytic near the origin and the singu- 
larity nearest to the origin is a pole of the order one. 
The pole f, and its residue 6_, are found to be 
| 0 
i. 2 sin?(ty) 


(23) 


b_, w~'2-de-%sin (hy). (24) 

The parameter y¥ in (22) depends on the direction 
of the incident wave and on the direction of the point 
of observation such, that y (and consequently fy) 
tends to zero when one of the boundaries of geo- 
metric optics is approached. The asymptotic ex- 
pansion of the integral (21) for large yp can then be 
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obtained by means of (18), (19), and (20): 


Tay) = >Se,T(n+4)(yp)-"~? 


n=0 


—he"""""GY) Erfe[(2yp)? sin y . (25) 
The coefficients ¢, are determined by (10) with 


g*(t) g(t) + w7'272e- sin (Ay) (t+-2 sin®y/2)~! (26) 
and g(t) is given by (22). The function g*(¢) fulfills 
the conditions (a) and (b) in the region —r<argt< x. 
The asymptotic expansion (25) holds therefore in 
—3n/2<arg (yp)<32/2 and is valid also for y=0 
(in which case all the ¢, vanish) and assumes in this 
case the correct value of (21) 


I (a,0) -e-o, (27) 


The results produced here are equivalent to those 
in [6,7]. 
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Principal Submatrices of a Full-Rowed Non-Negative 
Matrix’ 


K. Goldberg 


April 2, 1959 


It is shown that one of the principal submatrices of an incidence matrix without zero 


roWs is & permutation matrix. 
whether a non-negative matrix is nilpotent. 


We begin with a necessary flock of definitions. 
All matrices in this paper are square matrices. The 
principal submatrices of a matrix are the matrix itself 
and those submatrices obtained from it by repeatedly 
striking out a row and the column of the same index. 
The leading principal submatrices are those obtained 
by striking out exactly one row and its corresponding 
column. The principal minors are the determinants 
of the principal submatrices. 

A zero vector is one with every element equal to 0. 
A full-rowed matrix is one with no zero rows. A non- 
negative matrix is one whose elements are all non- 
negative real numbers. An incidence matrix is one 
whose elements are either 0 or 1.) A permutation 
matrix is an incidence matrix with exactly one 1 in 
each row and each column. If some power of a 
matrix vanishes the matrix is called nilpotent. We 
shall prove: 

THroremM. At least one of the principal submatrices 
of a full-rowed incidence matrix is a permutation matriz, 

Corotiary 1. Jf N is a non-negative matrix with 
a full-rowed principal submatrix of order m, then for 
some positive rm the trace of N’ is at least rn’, where 
n is the least positive element of N. 


Coro._uary 2. Each principal submatriz of a nilpo- 
tent non-negative matrix is nilpotent and has both zero 
rows and zero columns. 


Both the theorem and its first corollary hold with 
“row” replaced by “column”. The theorem relates 
to similar results on systems of distinct representa- 
tives.* 

The main effect of corollary 2 is to provide a simple 
algorithm for deciding whether or not a non-negative 
matrix is nilpotent. If it has no zero rows it 1s not. 


1 The preparation of this paper was supported (in part) by the Office of Naval 
Research 

27H. B. Mann and HJ 
Monthly 60, 397 (1953 


Ryser, System of distinct representatives, Am, Math. 


A 
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From this an inspection method is described for determining 


pplication is made to graph theory. 


Otherwise the zero rows and the columns of the same 
index are crossed out and the process repeated on the 
smaller matrix. If none of the matrices examined in 
this process are full-rowed, then all the principal 
minors of the matrix are 0. The characteristic 
equation is thus of the form «"=0,° and the matrix is 
nilpotent. R.Marimont* first described this process 
for incidence matrices. 

That the principal submatrices of a nilpotent non- 
negative matrix are also nilpotent is known. We 
have included this statement in corollary 2 because 
a new proof follows immediately from the remaining 
result in this corollary. 

We prove the theorem by induction. It is cer- 
tainly true fora matrix of order 1. Assume that it is 
true for matrices of order n-1. 

Let A=(a,,;) be a full-rowed incidence matrix of 
order n. If any of its leading principal submatrices 
is full-rowed then we are finished because “principal 
submatrix” is a hereditary property. Therefore, 
we assume that each of the leading principal sub- 
matrices of A has at least one zero row. We shall 
prove that then A is a permutation matrix. 

Let A, denote the leading principal submatrix 
of A obtained by striking out the row and column of 
index s. We may assume that there is an index 7, 
such that 
») 


0 > See a nm. (2.1) 


a))—0 for j 


Since no row of A is a zero row this implies that 


(2.2) 


ad, ,=1. 2.2 
If another leading principal submatrix A, also 
had as one of its zero rows a row corresponding to 
. . ° ‘rT ° 

the index 7, we would have a, ,=1 with t#s. This 


C.C. MacDurffee, The theory of matrices, p. 19 (Chelsea Pub. Co., New York, 
N.Y., 1946) 

‘R. Marimont, A new method of checking the consistency of precedence 
matrices, J. Assoc. Computing Machinery 6, 164 (1959). 








Therefore, to each s there 
which satisfies (2.1) and 
It follows that A is a permutation matrix, 
and the theorem is proved. 

For any non-negative matrix NV let V denote the 
incidence matrix of the same order as N, with 1 in 
every position that N has a positive element and 0 
in every position that N has 0. If 7 is the least 
positive element in N then N-»N is a non-negative 
matrix. Thus, if tr Y denotes the trace of LY, 


would contradict (2.1). 
corresponds a unique /, 
y) 2) 


tr N’>n’ tr N’ 


for every non-negative integer r. 

Furthermore, let N (n;;) be of order n, let 
S={1,2,..., n}, and let S, be any subset of S. 
Let A be the principal! submatrix of N with row 


and column indices in S). Then for any positive 
integer 7 we have 
» At. 
tr N’=5 ‘an, i, Misi - 2 MG, 
eS 
mi » r 
> Dy Ni, i, Miz igs © © M1 = tr M’. 
hee 


Now suppose .J7/ is a full-rowed principal submatrix 
of N of order m. Then JM is a full-rowed incidence 
matrix, and so has a principal submatrix 7? of order 
s < m which is a permutation matrix. Let 7 < S 
be any of the evele lengths of the permutation 
corresponding to ?. Then 


tr P’>r. 
Since P is a principal submatrix of 7, the incidence 


matrix corresponding to AZ, and AV is a principal 
submatrix of N, we have 
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trN>trM > qtr > 9 tr PF’ > ra’ 


with r < s < m. This proves corollary 1. 


Corollary 2 follows directly from corollary 1 and | 


the fact that a matrix is nilpotent (if and) only if the 
traces of its positive powers are all 0. 

Our theorem can also be stated as a result op 
directed graphs. A directed graph G is a set of 
points 2? PY: anda subset ??* Ce ee 
of the Cartesian product ? <P. The ordered pair 
(pP.qg) May be thought of as a directed line from p, 
called the starting point, to q 

We say that @ contains a simple loop if there is 4 
subset 1” =) pi.p), pr of P such that 
rr) 


P*qQ(P’ : (Pi.P2),(Ps.Ps), 


, ; 


(Pn —1sPn)s (PrsPi). 


That is to say, the points of 2?” have these lines 
between them, and no others. 

Our theorem is equivalent to the following state- 
ment: If every point of a directed graph ix a starting 
point then the graph contains a simple loop. 


- 


As so stated, this theorem is essentially one proved | 


by David Rosenblatt.’ 


5D. Rosenblatt, On the graphs 


tion matrices and stochasti¢e matrices, Naval Res 1957). 


Log. Quart. 4, 151 
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Zeros of Certain Polynomials 


A. J. Goldman 


(March 2, 1959 


Let P be «a real parameter. 


unit disk, if and only if 0< P< 2 sin x/(4n4 


1. Introduction. Let n be a fixed non-negative 
integer. The following problem will be solved: 

ProBLem. Find the set S, of all values of the real 
parameter P for which all roots of 


| P 
lie in the open unit disk \2\ <1. 


~nr 


0 (1) 


In view of the substantial literature relating the 
coefficients of a polynomial to the locations of its 
zeros, one would expect that some standard 
algorithm could be applied to so specific a problem. 
This seems to be only partly true; the appropriate 
algorithm is that given by the Schur-Cohen Criterion? 
which vields only an implicit characterization of S, 
as the solution-set of a system of n+ 1 polynomial 
inequalities in 77. It appears nontrivial to derive 
from this an explicit characterization of S,, and so 
we give instead an elementary self-contained solu- 
tion. The result is the 


P 


THrorem. S,, 0< P<2 sin r/(4n+-2) 
The problem arose in connection with the gen- 


erating function * 


ps ) 


Fr s) m9 
Pip ‘Ss 


1—s-+ ; 

of the recurrence times for runs of m successes in a 
sequence of Bernoulli trials with “success probabil- 
ity’ p. Rigorous justification of the usual probability 
theory manipulations of the power series for F(s) is 
easy if F(x) has no singularities for s <1, and this is 


In fact true and is equivalent (upon setting z=1/s) to 
the assertion that all roots of 2"*'—2"+ (1—p)p"=0 
lie in 1 for 0 p<1. The last statement shows 


that S, ‘neludes the interval {P| 0< P< n"/(n+1)""*")5, 
and it Was natural to inquire whether this expression 
gave S, exactly. It follows from our theorem that 
S, is larger than this interval for all n>0, and is 
approximately 'sme 4.27 times as long for large n. 


Marden, The geometry of the zeros of a polynomial in a complex variable, 
Am. Math. Sec, Math, Survey No. 3 (New York, N.Y., 1949). 

Ibid, p. 152. 

Feller, Probability theory and its applications, p. 265 
Inc. , New York, N.Y., 1950 


John Wiley & Sons, 


It is proved that all roots of 2"*! 
2). 


f' (0) 


"+ P=0 lie in the open 


2. Solution. In this section x, r, and @ denote real 
variables obeving r>0, 0<@<2r. Since the theo- 
rem stated above is obviously true for n=0, it is 
assumed that n>0 in what follows. It is convenient 
to define 


f(0)=sin"né sin 6/sin"*' (n+ 1)8, 
A r"—az"""| |g >1}, 
B=(0 0<60< 9, sin nO>sin (n+1)0>0}, 
C=(/(0) @inB 


We can omit the easy proof of 


Lemma 1. A=/{P| P<Oor P>2 if n is even, 
A P| P<o if n is odd, 
Lemma 2. C=:{P| P>2 sin x/(4n+2) 


Proor. (a) From the formula 


2nsin@sin n&(1—cos(n+ 1)8)] 


1)@ 


sin né) 4 
sin"*? (n-+ 


sin"! né[(nsin @ 


we conclude that /(@) is increasing on each sub- 
interval of B. 

(b) Suppose sin(n+1)@*=0 at a left endpoint 6* 
of some maximal subinterval of B. If sin né* #0, 
then sin 6*40 and so {(@*+)=(+ @), contradicting 
(a). If sin né*=0, then @*=0, contradicting the 
requirement that 


sin n(@* +6) >sin(n+ 1)(@*+-4) 


for all sufficiently small 6>0. Thus the supposition 


is untenable. 
(c) We next apply the identity 


sin né—sin(n+ 1) 2 sin 56 cos 5(2n +1) 
to obtain 
, e : 
B=(0 0<0<r, sin(n+1)0>0, coss(2n+1)0<S0}. 





Consider now any left endpoint @* of a maximal | 37/(2n+-1), but sin(n+1)@ changes sign earlier, 


subinterval of B. According to (b), we must have | g/(n+1). Thus @=2/(n +1), and so f(@—)=(+e)? 
This fact, together with (a) and the results of (¢ 
sin (né*)=sin(n + 1)6*, (2) completes the proof. 
so that — , . Our final lemma gives the motivation for lemma: 
F(0*) =sind*/sin(n + 1)0*. (3) | 1 and 2; the three lemmas together immediate) 
ne a ; imply the theorem stated in the introduction. i 
In fact, the conditions determining such an endpoint 
are, in addition to (2) and 0<@*<, that LemMMa 3. S, os the complement of AUC. 
sin(n+1)0* >0 and sin n(@* +6) >sin(n + 1)(@* +6), Proor. (a) Clearly eq (1) has a real root outside! 
a ae | the disk 2 <1 if and only if ? is in AL. 
for all sufficiently small 6 0. Equivalently, @* is (b) Next we observe that s=r exp(/@) is a nonregl! 
such an endpoint if and only if root of eq (1) if and only if 
l " SI mt! oi , la 
cos 5 On +-1)0*—0. (4) r” sin n@—7 in(m—- | (), (9 
sin 640, (10° 
me — (Pn « * 1 5). ree -_ . : 
~—" 2 si li ISM (6 as above) (9) r” cos n6—r"*! cos(n+1)0=P. (11 
. . ' 
sin(n+1)0* >0. (6) Now (9) and (10) are equivalent to 
' 
0<ot<x. (7) sin né/sin(n+1)0=r, (12 
; , , = ‘ ee a 5 
The points obeying (4), (5), and (7) are precisely and (11) and (12) are equivalent to : 
the points por 
f(ay—P. (13), tabl 
6;=(4j)+-1)x/(2n+1) (O<2j;<2n); we ; a 
Thus eq (1) has a nonreal root outside the disk |2)<j , prot 
these points also satisfy (6), since ~ if and only if 7 lies in the set lat 
. is Si 
sin(n + 1)0;=cos(4) + 1)r/(4n +2) f(8) sin nO/sin (n+1) 621), Mo 
(47+1<4n+2), (8) jae . . } not 
which (by considering the change 6—>27—9@) is readily tl 
' : : . Z aoe eS a é. ee wil 
and so the @,’s are precisely the left endpoints of the | Seen to be identical with the set ; 
maximal subintervals of B. From (3) and (8) we . = 
have f(0) sin nA@>sin (n+1)A>0). a 
Pr 
— , Base, 
{(0;)=2 sin(4) L1)r/(4n+2): (¢) By (a) and (b), S, is the complement ot W 
‘ } =), . NBs 
J 
together with (a), this shows that /(@) reaches its AU‘ f(8)) sin nd>sin (n +1) 6>0 ; 
minimum on PB at @=@—7/(2n+1), this minimum ; on fi. cl ; 
being By lemma 1, this union is identical with 


min ('=2 sin r/(4n+2). ; ; 
AU f(@) sin n6>sin (n+1) @>0, f(@)>0 


(d) Finally, we seek the right endpoint @ of that 


maximal subinterval of B of which 4 is the left end- which (by the form of f(9)) is just AUC. 


; l 
voint. After %, cos 5 (2n+1)@ first changes sign at , : 
, 2 . Wasuineton, D.C. (Paper 63B1-4) 
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Tables of Transport Integrals: A Supplement 


William M. Rogers,’ William J. Hall, and Robert L. Powell 


April 23, 1959 


were previously published in NBS Circular 595 for values of n from 2 through 17. 


e2"dz 
1) 


(e 


In this 


paper, values are given for the integrals where n is 18 and 20 and where zs ranges from 0.2 to 


values were calculated on the IBM 650 computer 


Tables of values of the transport integrals, defined by 
J, (r) 
50.0 in steps of 0 2 
Soon after the publication of the Tables of Trans- 
port Integrals? an unanticipated need arose for 


tables of the integrals for n equal to 18 and 20. The 
integrals with even values of » are useful in the 
problem of inverting heat capacity results to obtain 
lattice frequency spectra for solids. The technique 
is similar to the method of moments developed by 
Montroll.! The integrals with odd values of n are 
not used in this method, and since there was no other 
need for the integrals with » above 17, the integrals 
The additional 


with n=19 were not caleulated. 

Present address: Headquarters, Strategic Air Command, Offutt Air Force 
Base, Neb 

William M. Rogers and Robert L. Powell, Tables of transport integrals, 
NBS Cire. 595 

J. A. Morrison, private communication 

E. W. Montroll, J. Chem. Phys. 11, 481 (1943 
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using the already existent programs and auxiliary 
numbers. The computed values of the integrals are 
presented in the following tables. 

The mathematical formulation, analyses of errors, 
and tables of auxiliary numbers were discussed in the 
previous publication. For « below 10, interpolation 
should be performed by the technique of using an 
interpolating polynomial in a table of divided dif- 
ferences of the logarithms of both the arguments and 
the integral values. For + between 0 and 4, a fifth- 
degree polynomial should be used; between 4 and 10, 
a third-degree polynomial. For z above 10, straight- 
forward Lagrangian interpolation should be per- 
formed. For + between 10 and 20, a 4-point formula 
should be used; between 20 and 35, a 3-point; and 


above 35, a 2-point. 


(Tables follow) 





Rez 7N 


TRANSPORT INTEGRALS: Jp(x) = | =, dz; n=18 ,20 
































(e*-1)? 

x 1/x Jig (x) Joo(x ) 
COEF. EXP, COEF. EXP. 
e2 5.0000 7¢68716 | -14 2e75110 | -15 
oe 22e5V000 9e9861ll1 i- 9 12642940 | - 9 
06 1¢6667 9269389 |- 6 3012152 | - 6 
8 1262500 1e¢26317 | = 3 7022945 | - & 
120 1290000 5046290 | = 2 4288381 | - 2 
le2 e 8333 1217362 ) 1.51031 Q 
le4 e 7143 1255318 1 2e 71939 i 
16 06250 12643986 2 3229112 2 
128 2 5556 1¢01589 3 2093725 3 
220 e 5000 5077334 3 2205955 4 
2e2 04545 20e75197 a 1218710 5 
204 04167 1613372 5 5281587 5 
206 e 3846 4212961 5 20248433 6 
208 03571 1¢35385 6 9243812 6 
320 0 3333 4005157 6 3023956 7 
302 03125 1611950 7 1201753 8 
304 0294) 2288300 7 2295532 8 
306 2778 60e97382 7 8e00642 8 
328 e 2632 1259495 8 2203807 9 
420 e 2500 3046812 8 4.90503 9 
4e2 02381 7e¢20417 8 1212206 10 
4e4% 02273 1¢43552 9 2245095 10 
4e6 02174 2¢75377 9 5013255 10 
408 02983 5210162 9 1.03403 11 
520 22000 9215275 9 2201034 ll 
52 21923 1259414 10 3.78208 ll 
504% e 1852 2¢70137 10 6290194 il 
506 01786 4046243 10 1622443 12 
508 01724 7¢19868 10 2211576 12 
6eV 01667 1213584 ll 3056723 12 
6e2 e 1613 Le 75544 ll 5087786 l2 
60% 01563 2266088 ll 9047887 id 
6e6 01515 3096048 ll 1.49801 13 
6e8 e147] 50 19463 ll 2032279 13 
720 01429 8e 34235 ll 3453773 13 
Te2 e 1389 1218285 12 5229775 13 
704 01351 1265316 12 7280755 13 
76 21316 2e27919 12 1.13335 14 
7¢8 01282 3210193 12 1.62174 14 
829 01250 4017019 12 2228921 14 
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} 


) * e772" ‘ 
TRANSPORT INTEGRALS: Jpy(x) = / (e?-1)3 dz; n=18, 20 (contd) 
(e) 











} e7 -|)? 
x 1/x Jig (x) J2o(x) 
COEF. EXP. COEF. EXP, 
| 8e2 201220 5054145 | 12 3018984 | 14 
84 21190 7028248 | 12 4.39041 14 
806 01163 9047008 12 5697241 14 
8o8 01136 1e21916 | 13 8203411 14 
940 ellil 1255454 | 13 1.06928 | 15 
9e2 21987 1296411 13 1.40870 15 
904 2 1064 2245996 | 13 1.83786 | 15 
906 0 1042 3005533 | 13 2237590 | 15 
928 e 1020 30 7644] 13 3204320 15 
10,0 2 1000 4060253 13 3286510] 15 
10.2 20980 5058606 | 13 4e86890 | 15 
| 104 209962 6673205 | 13 62008530| 15 
1006 00943 8 205832 13 72054830] 15 
' 1068 00926 92058320! 13 9.29500] 15 
1120 eV9U9 1¢13254 14 1213656 16 
lle2 20893 1633040 | 14 1238045 16 
lle4 ©0877 1655380 14 1e66582 16 
lleé 20862 1080461 14 1099760 | 16 
1168 00847 2208468 14 2238115 16 
12.0 20833 2039576 | 14 2.82181 16 
1202 200820 2273954 | 14 3032528 16 
1204 20806 3011761 14 3089742 16 
1206 00794 3053141 14 4054416 | 16 
128 20781 3098224 | 14 5027146 16 
1360 20769 4047120 14 6008533 16. 
1302 20758 4.99927 | 14 6.99174 | 16 
1304 00746 5056715 14 7299644 | 16 
1306 00735 6017539 | 14 9.10518 16 
1368 00725 6082427 | 14 1203233 | 17 
1440 00714 72051388 | 14 1216559 | 17 
1462 20704 8024402 14 1631077 | 17 
1404 20694 9201431 14 1.46831 17 
1406 20685 9682409 | 14 1.63858 17 
1468 09676 1206725 15 1.82193 17 
1520 20667 1215583 15 2201862] 17 
1502 00658 1e 24803 15 2022887 | 17 
1504 20649 1¢34370 | 15 2045283 | 17 
1506 2064) 1044264 | 15 2069057 | 17 
1508 00633 1054468 | 15 2094210 17 
1620 20625 1264960 15 3220737 | 17 
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x zZz5Nn 
TRANSPORT INTEGRALS: Jp(x) =f <*— dz; n=18,20 (cont'd) ’ T 














(0) (e7-1)? ) 
x 1/x Jig (x) J2o(x) 
COEF. | EXP. COEF. | EXP. 
L602 00617 1675718 | 15 3048623| 17 | 
16e4 eV610 1e¢86717 15 3e77849 17 
1606 20602 1097933 | 15 4-08386| 17 
1668 20595 22009340] 15 4.40200| 17 
17.0 20586 2e20911] 15 4.73249 | 17 
17e2 00581 2¢32619] 15 5007485] 17 
1704 00575 2044436] 15 5042855] 17 
1706 00568 2056336 | 15 5.79297 | 17 
1708 00562 2068289 | 15 6016747| 17 
1860 00556 2e80270] 15 6655134| i7 ' 
1802 00549 2092250] 15 6.94383 | 17 
1804 00543 3004204] 15 7034417| 17 
1806 00538 3016107] 15 7075153 | 17 | 
1868 03532 3027933 15 8.16507 17 
19,0 20526 3039658] 15 8.56393 | 17 
19e2 00521 3051262] 15 9200721! 17 
1904 00515 3062720 | 15 9643404 / 17 
1906 20510 3074015 | 15 9686352 17 
1908 00505 3085127 15 1.02948 18 
20.0 09500 3096039] 15 1.07269! 18 
2062 00495 4006735] 15 1011590! 18 
2004 00490 4017201 | 15 1615902| 18 
2066 00485 4027423 | 15 1020198! 18 
2008 0048] 4037390) 15 10624469 | 18 
2100 00476 4047093 | 15 1.628707/| 18 
21602 00472 456522 | 15 1.32905} 18 
2104 00467 4e65671 | 15 1237056! 18 
21e6 00463 4e 74534 15 1041152 18 
21e8 0459 4e831VU7 15 1645189 | 18 
2200 00455 4091386 | 15 1649159 | 18 
2202 20450 4099369 |} 15 1053058 | 18 
2204 0 0446 5007056 | 15 1056881! 18 
2206 00442 5014447 | 15 1.60622 | 18 
2208 00439 5021543} 15 1.64278 | 18 
2300 0 V435 5028346) 15 1067846 | 18 
2302 00431 5034860 | 15 1e71322| 16 
2304 00427 5041089 | 15 1674703 | 18 
2306 00424 5047037 | 15 1.77987 | 18 
2308 0 042U 5252709 | 15 1.81173 | 18 
2400 00417 5658111} 15 1.84259 | 18 
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d) ' TRANSPORT INTEGRALS: Jp( ef Si 


e727" 


dz; n=18,20 (cont'd) 














) e z_1)2 
x 1/x Jig (x) J2o(x) 
COEF. EXP. COEF. EXP. 
\ 24e2 094)]3 52063250 15 1287243 16 
24e4% 00410 5268132 15 1.90126 16 
I 2406 eV4U7 5072765 15 1.92907 18 
24ed 00403 5077155 15 1.95585 18 
2500 20400 5281312 15 1298162 18 
2502 00397 52065242 15 2290638 18 
2504 003946 5288954 15 2203014 186 
j 2506 e049] 5092455 15 2205290 18 
2508 09388 5095755 15 2207470 18 
' 260 eVJ3a5 5e¢ 98862 15 22099553 18 
2602 e4¥382 60017683 15 2211543 18 
2604 e379 62004526 15 2213441 18 
: 2606 0376 60e0710)] 15 2015249 18 
2608 e373 6009514 15 2216969 16 
27¢e9 eV37uU 6e11774 15 22168604 16 
2702 eV 368 6013889 15 2020157 18 
2704 0 365 6¢15865 15 2e21629 18 
2706 0% 362 6017710 15 2023025 18 
2708 00360 6619431} 15 2024345] 18 
2820 00357 60e21034 15 2025593 16 
28e2 200355 6022528 15 2026772 18 
2804 00352 6¢23917 15 2027885 18 
28e6 00350 6025208 15 2028933 18 
28e8 00347 6026407 15 2029921 16 
2940 00345 6e27519 15 2230849 18 
2902 00342 6028550 15 2e31722 16 
29 0% 29340 62029505 15 26 32542 18 
2906 00338 6¢30388 15 2¢33311 18 
2908 00336 6¢31205 15 2¢ 34031 16 
30,0 0333 6¢31960 15 2034706 18 
3002 eV33] 6032657 15 2¢ 35337 18 
3004 00329 6¢33299 15 2035927 18 
30065 00327 6033892 15 22¢ 36478 18 
3008 00325 6034437 15 22036992 ié 
31¢VU 29323 6034939 15 2037472 18 
3le2 eV321 6¢3540] 15 2237918 18 
3164 09318 6¢35825 15 2e 38333 18 
316 09316 6¢36214 15 2238720 18 
318 09314 6036572 15 2039079 18 
3260 00313 6¢ 36899 15 2039412 18 
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x ez 7" 


dz; n=18, 20 (cont'd) 











TRANSPORT INTEGRALS: Jn (x) = (e212 
x I/x Jig (x) J20(x) 

COEF. EXP, COEF. EXP, 
3202 e931] 6¢37199 15 2039722 18 
320% 09309 6e37474 15 2240008 18 
3206 09307 6¢37726 15 2040274 18 
328 209305 6037956 15 2240520 18 
3340 09303 6¢38166 15 2040747 18 
3302 e030] 6¢38358 15 2040957 18 
3304 e299 6¢38533 15 2e41151 18 
3366 04298 6¢38692 15 2e41330 18 
3368 00296 60 38838 15 2041496 18 
3420 00294 6638970 15 2e41648 18 
34e2 e%V292 6¢3909)1 15 2241788 18 
3404 e029] 6¢ 39200 15 2041917 18 
3468 e289 6¢393V0 15 2042036 18 
3408 eV287 6¢39391 15 2042145 18 
3520 e%286 6039473 15 2042245 18 
3502 03284 60639548 15 2042337 16 
3504 00282 6e39615 15 2e42421 18 
3546 eV281 60¢39677 15 2042499 18 
3508 eV279 6039732 15 2042570 18 
3620 200278 62e39783 15 2042634 18 
3602 0e%276 6¢39828 15 2042694 18 
3604 00275 6¢39869 15 2042748 18 
3666 00273 6639907 15 2242798 18 
3608 e272 6639940 15 2242843 18 
3700 00270 6¢39971 15 2042884 18 
37e2 00269 62e39998 15 2042922 18 
3704 009267 62¢400 23 i5 2042956 18 
3706 e266 6e40U45 15 2042988 i8 
378 0%265 6040065 | 15 2043016 18 
38.0 09263 6040083 | 15 2043042 18 
38e2 0%262 6040099 | 15 2243066 18 
38 «4% 02260 6040114, 15 2243087 18 
3806 09259 6040127 | 15 2043107 18 
3608 09258 6¢40139 15 2043124 16 
3940 eV256 6040150 15 2¢43140 18 
39e2 00255 6040159 15 2043155 18 
39 «4 00254 6¢40168 15 2243168 18 
39«6 09253 6040175 15 2243180 18 
398 e025] 640182 15 2243191 18 
40.0 09250 6e40188 15 2243201 ig 
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J) 











x ZN 
TRANSPORT INTEGRALS: Jn(xd ef SA, dz; n=18, 20 (cont'd) 
A ms 
x 1\/x Jig (x) J20(x) 

COEF. EXP, COEF. EXP. 
402 03249 6240194 15 2e43210 18 
4064 04248 6240199 15 2243218 18 
4026 00246 6040203 15 2043225 18 
4028 0245 6e40207 15 2043232 18 
41.0 09244 6e40211 15 2043238 18 
41le2 09243 6040214 15 2043243 18 
414 09242 6040217 15 2043248 18 
416 04240 6240219 15 2043252 18 
41.8 e239 6¢40222 15 20 43256 18 
4220 03238 6240224 15 2043260 18 
4202 00237 60e40225 15 2043263 18 
426% e09236 6040227 15 2043266 18 
4226 09235 6240228 15 2243268 18 
4248 09234 6040230 15 2043270 18 
4329 eV233 6240231 15 2043273 18 
4302 0923) 6040232 15 2043274 18 
430% eV23U 60¢40233 15 2243276 18 
4346 eV229 6240234 15 2243278 18 
4368 09228 6040234 15 2243279 18 
4440 eV227 6040235 15 2243280 18 
44.2 04226 6¢40235 15 2243281 18 
44064 23225 6040236 15 2043262 18 
4406 0oV224 6040236 35 2e43283 18 
4468 009223 6240237 15 2043284 18 
4560 eV222 6040237 15 2043285 18 
4502 09221 6¢40237 15 2043285 18 
4504 09220 6¢40238 15 2043286 18 
45066 eV219 6240238 15 2243286 18 
4568 00218 6240238 15 2043287 18 
4620 00217 6e40238 35 2043287 18 
4662 e0216 6040239 15 2043288 18 
4604 eV216 6040239 15 2043288 18 
4666 eV215 6040239 15 2243288 18 
46e3 eJ214 60e40239 15 2043288 18 
476% 0%213 6040239 15 2243289 i8 
472 e212 640239 15 2243289 18 
474% eY%211l 620240239 15 2043289 18 
4706 29210 6240239 15 2243289 18 
4728 209209 62040239 15 2043289 18 
4BeV 09208 6040239 15 2043289 18 
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dz ;-n=18,20 (cont'g) 











TRANSPORT INTEGRALS: Jn(x) =f <— 
o(e —1) 
x l/x Jig (x) J2o(x ) 

COEF. EXP, COEF EXP, 
482 29207 6040239 15 2043290 18 
4804 00207 60e40239 i5 2043290 16 
4826 00206 6240240 15 2043290 18 
488 09205 6040240 15 2043290 16 
49-20 00204 6040240 15 2043290 18 
4902 eV203 6240240 15 2243290 18 
4904 09202 6¢40240 15 2043290 18 
4926 00202 6040240 15 2243290 18 
4928 e020) 6040240 15 2043290 18 
50.29 eV200 6040240 15 2243290 18 
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(Paper 63B1-5 


Journal of Research of the National Bureau of Standards—B. Mathematics and Mathematical Physics 


Vol. 63B, No. 1, July-September 1959 


Lens Design: A New Approach’ 


Orestes N. Stavroudis 


(February 18, 1959 


This paper describes a new method of defining the total aberrations of an optical system 
and its application to lens design. A system of partial differential equations defining optical 
image formation is written in the form of 4 by 4 matrices and is applied to the derivation 
of lens aberrations. It is shown that the form of the matrices facilitates the computation 
of these aberrations and that equations of condition on the lens parameters can be derived. 
The use of this method is illustrated through its spplication to a simple, well-known problem. 


1. Introduction 


In the course of designing a lens, one invariably computes aberration functions whose 
magnitudes reflect: qualitatively and quantitatively the performance of an optical system. 
Unfortunately these quantities are rather difficult to calculate, usually requiring a long program 
of ray tracing, and, more importantly, are not readily correlated with the lens parameters. 
Usually a lens designer must rely on his experience when determing whether a change in 
a certain curvature, thickness, or index of refraction will affect a given aberration. 

In this paper a method of attack will be described which will at least mitigate these dif- 
ficulties. The application of this method results in the definition of a set of aberrations to- 
gether with a systematic scheme for their computation expressed as equations in which the 
system parameters appear. Yet the method ts sufficiently flexible to meet the special require- 
ments of any lens. For example, one system of aberrations can be defined for use with a copy 
lens and a second system can be constructed for application to an afocal system. 

Consider an optical system consisting of a system of lenses and fixed reference planes in 
object and image space which are not necessarily object and image planes. Right-handed 
coordinate svstems are chosen in object and image space, the coordinates of which are desig- 
nated av. y, 2, and w’, y’, 2’, respectively. We take the positive z and 2’ axes in the general 
direction of ray propagation and v-y plane and 2’-y’ plane coincident with the object and 
image reference planes. Let the reduced direction cosines of a ray in object and image space 
be & 9. Gand &, n’, ¢’, where &4-*n9+ CaN? EP 4+? + FP = N”, N and N’ being the indices of 
refraction of object and image space 

Associated with each ray passing through the optical system are two sets of four independ- 
ent coordinates, , y, & , and 2’, y’, &, 7’; for a given optical system either set completely 
determines a ray. [nan obvious way a given optical system can be said to define 2’, y’, &’, and 
n’ as functions of x, y, & and». In what follows we tacitly assume that in any optical system 
these four funetions are continuous and possess continuous first-partial derivatives. 

Suppose every ray emanating from an object pomt converges to a single pomt on the image 
plane. Then «’ and y’, considered as functions of the object space variables, are unaffected by 
a change of € andy. Thus «’ and y’ are independent of € and 7 and 


t,=2 Yi y,—0 (1) 
where the subscripts denote partial derivatives. The converse is also true; if there exists a pair 


of values for « and y for which (1) is satisfied, then all rays emanating from the point (a, y) in 
object space converge to a single point («’, y’) in image space. Thus the four equations of (1) 


This paper contains work done at the Imperial College of Science and Technology and included in a thesis submitted to the | niversity of 
London in partial fulfillment of the requirements for the Ph. D. degree The work was supported in part by the U.S. Air Force, 
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provide a necessary and sufficient condition for a point to be imaged perfectly. By a similar 
argument it can be seen that in a region where (1) is satisfied a necessary and sufficient condition 
for a lens to be free from distortion is that 


, , , 


t,=Y,=—™M, z. Ys Q, (2) 
where m is a constant and represents the magnification of the system. 

These two sets of equations provide conditions assuring perfect, distortion-free image 
formation. By the same token they may be taken as measures of the failure of an optical system 
What 


to achieve perfection; that is to say, they are measures of the aberrations of the system. 


follows is an elaboration of this basic idea. 


2. The Fundamental Optical Invariant 


The functional relationship described in the preceding section is a general transformation 
of which optical image formation is a special case. It is desirable therefore to impose restric- 
tions on the functions; to exclude from our attention those relationships that do not apply to 
optical systems. There are several ways in which this can be done. For the purposes of this 
study the most convenient method is through the use of Herzberger’s fundamental optical 
invariant. This leads to a system of partial differential equations which must be satisfied by 
the functions 2’, y’, &’, and 7’ if they are to describe the image formation of an optical system. 
Following Herzberger [1] these are written, 








Tety t+ Yr, —2,t2—Yy92=0 ) 
47S: Urn: I, YeNe I 
tq Urq 4 Se — Ys =O 

PE ee ae (+3) ’ 
LySe > YyNe — 2 Sy —YeMy =O 

tas . _ ae } 
LySn 7 Yun WySy YnNy I 

, / , , , Pa , , 
Le Sq Ye My — Ly Es — Yn Me (). 

It is particularly convenient to express these in matrix form as 
XS X J, (4) 


where 














and A represents the transpose of LY, 
It is important to note that X is a Jacobian matrix. An interesting corollary to (4) is that 


Ala, (5) 
or, in words, the Jacobian of the transformation representing optical image’formation is equal 


to unity. 
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It is interesting to note that this corollary is equivalent to a theorem proven by Straubel [2] 
as long ago as 1902. He applied this to several problems including photometry. More recently 
Toraldo di Francia [3] used this theorem in describing a generalization of geometric optics. 

The Y-matrix is particularly useful. First of all its elements consist of the partial deriva- 
tives of image space variables with respect to object space variables and therefore can be said to 
represent the geometrical aberrations of an optical system. Moreover, it appears as an indeter- 
minate in the equations derived from the fundamental optical invariant. Finally it can be de- 
composed into a product of .Y-matrices each representing a component of the lens. Thus, if we 
consider an optical system consisting of two lenses arranged so that the image plane of the 
first coincides with the object plane of the second then 


N=X2X, (6) 


where .Y, is the Y-matrix of the first lens, Y, of the second and Y that of the composite op- 
tical system. 

The proof of (6) follows at once from the properties of Jacobian matrices and the rules for 
differentiating a function of a function. Let the first lens be defined by the four functions 2, ¥, 
| E, and » of «, y, &, and 7 and let the second be defined by x’, y’, &’, and 7’ of 7, y, ¢, and ». The 
composite lens is then defined by 2’, y’, &’, and 9’ as functions of z, y, £,9. Then differentiating 




















, ed of 1 ry \=- 
r I-7-+Y- -7-£,+7-n,- 
| a z72 Y>Ur ySe re%: 
, i- ay, re ss { 
) —Let-Y- Yet Y-Eet », Cte. 
" y, Y>42 UU Yess Y>Nz, & 
ry’ 
Then 
P . i Be , , ‘ = ad - > 
i 2 Ze a 
(a: Vy dz dy § . ss “S Ty Ly Yt Ly 
, , y. , , , —_ = —_ om 
‘a ~ < /- 
Ve Uy Ye Un mY YE Y, Yr Uy Yt Yn 
’ X 
, ’ 9 » t’ ¢f tog - - - — 
- gt ot » ? 2! pf pt i ee 
Sz Su Se Sn Sz Sy SE Sy Se Sv & &, 
se | ee a a ae 
union, me md 1 2 UD ty Me Med 
XeA1- 


Clearly this can be extended to lenses consisting of any number of components. 


3. Transformation of the X-Matrix 


So far we have been dealing with the variables x, y, &, and ¢. It is convenient to apply 
a transformation, T, mapping these variables into a new system of coordinates, say a, 8, y, and 
6, given formally by writing the latter as continuously differentiable functions of the former. 
The same transformation 1s applied to both the object and image space variables. 
If we designate the J/-matrix by 
@), & a), a; 


B. Bs B, B; 
M= sibs hig ? 
Ya Ys 77 Ys 


, 


bo 54 6, Os 


then 
M=A’'XA-"', 
where 
Qty Gy Ar Ay» 
no B, B, Be B, 
Yr Yu Vt Yn 
L5, 6, 5 by 
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A’ represents the same matrix but in image-space variables. Applying this transformation 
to equation (4) we have 


(A-'"XA’) A’ JA’-"(A'XA-)=A-' JA 
which reduces to 


M(A’-'JA'-")M=(A’-!JA-'). (8) 
If the transformation, T, is restricted by the equations 


AJA=J, A'JA’=d. (9) 
then (8) becomes 


MJ M=d. (10) 


which is identical in form to (4) except that the variables are changed from s, y, € and » to 


a, 8, y, and 6. Thus we can change our variables in any way that is convenient, provided 
that (9) is satisfied, without adding to the complexity of the fundamental optical invariant. 


4. Rotational Symmetric Systems 


In the case of rotationally symmetric systems we are in a position to make a particularly 
advantageous choice of one of the new variables. In such a system the skewness invariant 
holds and 

p=©r 7 yég rn— yé. 


e. 2 


By letting Y=rn— yé=p, and y’=—2'n’ —y’t’, the M-matrix degenerates to 


M 


When this is substituted into (10) we obtain the transformed fundamental optical invariant 


Os 8,64 8,6, (11) 


8. 6; B; é, 0) ( 12) 


We may consider the first two equations of (12) as a linear, homogeneous, simultaneous 
algebraic system in 8, and 6,. A necessary and sufficient condition for a non-trivial solution 
to exist is that the determinant of coefficients must vanish identically. However, the third 
equation of (12) shows that this determinant is identical to a non-zero constant. We conclude 
that B2—6,—0. Therefore the first equation of (11) reduces to aj=1. The remaining equa- 
tions of the transformed fundamental optical invariant become 


a3 = Bydg— B36, 
, 


Set Pea 
Os B10; B56», 
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and 


856; —B;6,—=1. (14) 


Solving (13) for 8, and 6,, using (14), we obtain equivalent expressions 


, —. > 
By = B30 — Brag 
(15) 


> 


, = a a 
p 6.05 0; Ag. 


Thus, under the assumption of rotational symmetry the \/-matrix reduces to 


H Qs a’ a; ) 


p 


oh * B; Bi, B; 





M ; . . (16) 











5. Further Simplification 


Suppose we have a rotationally-svmmetric optical system defined by the functions a’, 
3’ and 6’ of a, 8, p and 6, consisting of two components, the first defined by a’, 8’, and 6’ of 


a. 8, P, and 6 and the second by @, 8, and 6 of a, 8, p and 6. Of course p and P are identical 
numerically but play different logical roles which must be kept distinct. As indicated in (6), 


l a; a, a; | ae ae a: ] Qs ay Qs 
0 8; 8, B] 10 Bo B B11 0 Bs B, Bs 

7 oe . (17 
oOo i 9 © it @ 0 0 1 0 


65 


~ 
- 
a) 
~ 
So 
mo 
U 
I! 


p 


By interchanging the two right columns and the two bottom rows of each matrix in (17) we 


obtain 
l a, a of l at a! al) (1 ae Gs &, 
fe 6 p 
0 B; B B, 0 BB B11 0 Bs Bs B, 
0 65 4; : 
0 0 


sf i ey ->cfy (18) 
6; 46, 0 0- 0; 0- QO 62 3b; 6, 
Oo i 9 86 o> I oS 0 @ I 
Now let 
A’ (a3 ,a5), BR’ (a), 
c7 (7), pr=(*), 
636; b, 


with A, B,C, D and A, B,C, DP defined similarly. Then, by partitioning the matrices of 
(18) we obtam 
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¢ A ly’ r Cc D ¢ (' dD (19) 


0 ‘ 
where e=(5) and f=(0,0). Then from (19) we can obtain 


C’=CC 
or 
Bs 8; 83 B5\ (Bs Bs (20) 
, , , a? - 4 » 
63 65 55 05 63 65 
It has been shown that a consequence of (11) is that a’,=1, whence 
, , : 9) 
a’ =a-+X' (p, p, 6). (21) 
Applying this relation to the two components of the svstem under consideration, 
a’=at+NXB, P, 4) 
a =a+X(B, p, 5). 
Thus 
d (22) 


a’ —a=d’=)d+ X. 
— : : ; , ; ; 
Thus the geometrical properties of a rotationally -symmetric optical system can be described 
completely by the system of equations 


" 


a’ —a=>) X, (23) 
. i=] 
and 
8 8; n [{ O8,/OB;-, O06 ;/06;—; 
—= II , (24) 
53 6; i=] 06 OB, l O06, 06, l 
where A ;., is a function of 8;, p and 6;, where B’=8,, 6’=6,, and where B= By, 66). This 


is a system of five scalar equations of which four are independent by (14). 

If we replace the left members of these equations by functions representing a perfect 
optical system we are left with equations of condition on the parameter occurring in each of the 
components of the system. Alternatively a measure of the failure of the equations to be 
satisfied can be taken as a measure of the aberrations of the optical system. 


6. Functions for a Perfect Optical System. 


A possible choice for the variables, suitable for work with an optical system with object 
and image planes located a finite distance from the lens is 


a=aretan y/r, ) 





p=yi'+y’, > > 
(2-)) 
TET YN 
0 = =* 
yie+¥? J 


It can be verified that these satisfy (9) when y= p. 


By integrating the equations for a perfect optical system (1) and (2), we obtain 2’ =mze, 
y’=my, the constant of integration being eliminated by a suitable choice of coordinate axes. 
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From this 
a’=a, and B’= m8, 


whence 8’s=m and 8’;=0. Substituting this in (14) we obtain 


mo’; a 
Integrating: 
5’ =5/m+y, 

where ¥ is an arbitrary function of 8. Thus for a perfect optical system described by the 
variables in (25) we have 

a’ —a Dr, 0, (26) 

i 

and 


(27) 


i 08,/d8, 08,/26,—1° 
(, 1/m = (a5 a8, 05, a) 


as our equations of condition. It can be shown that ¥s when 8=0 equals the reciprocal of the 


1 


focal length of the perfect lens. 
The complexity of the problem may be reduced somewhat by the following considerations. 
Suppose (27) were satisfied for meridian rays only. Then 


B’ = mB-+-«,(B,6,p) 
5’ =6/m +-e,( B,6,p) 


where ¢=e@—0 when p=0._ If (26) is subsequently corrected for skew rays then a’s=a’s=0, 
which by (15) implies that 8’,—6’,—0. Thus the derivatives of « and & with respect to p 
are zero. Therefore these two functions are independent of p. Since «, and & vanish for a 
particular value of p, they must equal zero for all values and (27) 1s satisfied for skew rays as 
well as meridian rays. Thus, the solution of (27) for meridian rays and the solution of (26) 
for skew rays assures the correction of the optical system for skew rays. 


7. Transfer and Refraction Calculations 


The ultimate components of an optical system will be considered here to consist of (a) a 
transfer across a space between two planes and (b) a transfer from a plane (which will be called 
the reference plane) to a refracting surface, refraction across the surface and a transfer back 
to the reference plane. The ray tracing formulas for the transfer are 


Y y-+tn/§ 
oad t 
a s 
n =n 


where ¢ is the distance between the two planes, whence 
tp 
=a-+ are tan —- =. 
a a B+ 16 
a+ are tan p/pé—are tan p/p'd’, 
(28) 
a’ =, (B+ 18/¢)?+ (tp/Be)*» 


5/ 6 (B+ 16/¢) +-(tp/BS)p/B 
Lf . 
y (B+- 16/6)" + (tp/Bo)* , 
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For meridian rays the latter two formulas are 


p’ —(B-+t6/¢) ) 
} 29) 
J 


where @=N°—6&. Accordingly the matrix representing transfer between two planes separated 
by a distance ¢, in a medium with index of refraction WV, is 


l N2t/2° 
i T ° (30 
0 ] j 


The ray tracing formulas for a transfer from a plane to a refracting surface, refraction 
across the surface, and a transfer from the surface back to the plane are, respectively 


roar+ 2E/i ’ 

7] Yr =ni 

t’—t+ro 

n’=n+ yo { , 





where = is a function of 8 denoting the refracting surface and where 


a=-a-+ arctan ——= 


B= (B+28/¢)?+(Zp/Be)? > 








@  5(8+26/6)+(p/B)(=p/Bo) 
5 = = 
y (8- 20/6)° + (zp/B¢)* J 
a ‘ 
0 o> Do 
+ 4 
\ ¢ fete) 3 3 
: > 
, == -p ~ ) ) 
a a—arctan — —_ =a—arctan = aretan 
B°¢’ — =B6 Do 36 
3’ =-y (B— 28’ /¢’)? + (zp/Be") 4 “ 
8 (B23 /t") —(p/B) (Spl Be’ 
0 — — — 
Vip = 6" /% *-+. (2 p/Be’)* J 


' 


Thus for skew rays 


, p Pp P 
a a—arctan = >-+arctan ———arctan 5=+arctan 
po 86’ Bd 36 


ry pB(d—6") | 
a—arctan —>j-Tarctan =~=s arctan 
po B*66' +p 
) I@D . ; an 
a—aretan I _ arctan —— Pa = arctan Ss ov 
B 0 bdo 4 6- P, B-) Bé 


For meridian rays the formulas for 8 and 6 reduce from those given In (33), (34), and (35 


to 
p= B+ 26 ¢ ; 
4 
0 6 ; 
p’=, 
‘a ae Ss 
0 6 Do, 
] 
nha 


where 4 NV? o- ana c’? N*’ 6” ‘ 


By differentiating the first equation of (37) with respect to Band 6, recalling that = isa function 
of 8, we obtain 


B3= 1 + 23936/¢ 
re) 2 ( 238 6 4 =6 4 
so that 
33= C/A 
3 C/A) \ 2/e 
where 
A= <9 


Thus 


636 


V differentiating the equations of (38), 


B’3 8B’; | 0 
Mi» _ . “[-  - = ; - = z= i 
6363 bc > 3a? 3A 1 +-Bd/6-— TBC) A’¢ 
Where = ¢’—776’. From the definitions of A and A’ and from (34) it can be shown that 


| + BO(6—> 276) / A’ = Ae’) A’e, resulting in a simplification of m,. 


HUSUSH§—IU 4 39 








Again, differentiating (39) we obtain 


When these three matrices are multiplied together, we have 


MR=MeMpMe 


1 EN (aiiar | 0 PEN Te! (40) 
0 | BES’2 5530/2544’ Ag’/A’t 0 | 


Note that the first and last factors are formally identical to mz and can therefore be combined 
conveniently with an adjoining transfer matrix. 
When the refracting surface is a sphere centered on the reference plane 


2 Wi 2B. 


7 ( 


The minus sign assures that the sphere is concave to the right when ¢ is positive, in accordance 
with the usual convention. Then 





I —2N'/¢"\ ( A’ ¢/At’ 0 | =N2/¢ 
ma=(g M co'o/AZAA’ —s At’ vel 0 | ) (41) 
where 
ee > 
A= (66+ 2¢) 
as = | so 
A’ =(p5’ + 2¢’) y* = (N?—N"*)> p (42) 
C"< 
o c?z(A—A’). i 


8. A Simple Illustration 


As an example we consider the properties of a single refracting spherical surface. Con- 
sider an optical system consisting of an object plane, a spherical surface separating media of 
refractive indices NV and N’ and an image plane. Let the object and image planes be located 
distances ¢ and ¢’ from the center of the spherical surface. Figure 1 shows the case where ¢, 





Fieurre | 
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t’, and ¢ are positive. From (27), (30), and (41) the equation of condition for meridian rays is 
‘pk 0 ) 
V3 1/m 


| Nt’ /¢?\/1 —2ZzN"/¢~" A’¢/Ag’ ZN7/8\/1 N?t/¢' 
at l Yo l “ Meerpjetstan’ Pe i: l BF l ) 


For skew rays, by (26), (28), and (36), 


a’ —a=—arctan P +arctan P _aretan = po as =n), (44) 
B'6 B6 Bio+ (0° +-p*/8°) 


Expressing (43) in scalar form we obtain four equations. 
Ac/dc’ +N(t’ —2) to/e2"AA’t"@=m_ 

tool 4 | JE QA AP 
Fe pT AL = Ye 


N?(t+2Zz)d’ , N’(t’—z)A , N*2(t+2)N"(t’—2)@ 


~ Te TT =0 
APE A’gg” 72° Ad’ ¢?¢" 





At’/A’t + N2(t+2) ¢'o/222A4' =1/m. J 
We will pay particular attention to the third of these, which we write in the form 


N2(t+ 2)A"¢’ + N’2(t’ — 2) A2e+ (1/2?) (t+ 2)N"(t’ —z)o=0. (46) 


o’ = A’z= BS’ +20 =y oe? —(1/e)(N?—N”). 
Also, 
= +29, 
where @ e(a¢—a’) by (42). Substituting these vields 
eat {N’2(t’ —2)+ N2(t4+-2)] —N*%¢(t4+2)(N2?—N”) + N2(t4+-2) [Coz —N22+N "t’']o=0. 
This may be written as 


cote (Nt +N) + ¢(N2—N") [e202 — N2(t-+-2)] +N2(t+2)[e@et2 —N22+N"t’]o=0. 


Nt’ + N%t=0, (47) 
Then 
[¢(N2?—N’) + N2(t+2)¢][eo22 — N2(t+2)]=0. 
We thus obtain two possible cases represented by the two equations 
¢(N?—N”)+ N2(t+2)o=0 (48) 


eoe?=— N2(t4-2)=—0. (49) 








When 8=6—0, 


2 l/¢C.o N/e,o’ N’Je,c=N, and @=c 


“ 
a 


Substituting in (48) we obtain 
N(N?—N”) + N2(t—1/e)e(N—N’)=0, . 4 
N+N’+N(ct—1)=0, 


Applving this to (47), 
t’=N/N'e. 51) 


When B and 6 are not zero, it can be shown that (50) and (51) satisfy (48) for all values of 4 


provided that B—0. Therefore only the axial point on the object plane is imaged perfectly. 


This result provides the aplanatic points of a sphere [9] (see fig. 2). 
a 
} +| ~ 
r . T 
: tl 
a 
Picture 2 
On making the same substitutions into (49) we find that #—0, which when applied to (47 1 ow 
, , = = Wl 
vields t 0. As im the previous case, these are true for all values of 6 and 6 provided that we 
5=@. This illustrates the well-known fact that the center of a sphere is imaged perfectly on 
itself [4]. : 
—————_——— ch 
ry . . ° cn hl ° . mn 
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The integration of the storage differential equation at present is usually done mostly by 
| | ) y Dy) 


graphical or numerical procedures. 
tion is the subject of this paper. 


mathematically tractable expressions is introduced. 


is expressed in the form of a power function. 


storage y’ +cPi? —cy* 


shape and type of flow, and P being the inflow hydrograph. 
constant, and certain P 


0, P 
29)). 


equation for P 
to 3, eqs (12) to 
paper, 


1. Introduction 


The known simple balance equation, with the 
inflow minus the outflow in a time interval equal to 
the storage change in that interval is 


dS 


ie oe (1) 


which ts normally derived from the definition itself; 
with P the inflow, @ the outflow, S the storage 
volume, and ¢ the time. 

The two basic differential equations which govern 
the unsteady water movement through lakes and 
channels are the continuity and the dynamic (or 
momentum) equations: 


O(Ar) . dul 
“ex °) 
or or oh 
, as 3 
Stt ort I oe Ps \) 


where 

cross-sectional flow area, 

r—mean velocity of flow, 

z—length along the lake or the channel, 


‘time, 
g—earth acceleration, 
h—level of water surface in a cross sec- 
tion with reference to a fixed level, 
and 
S,— frictional slope. 
--shh 


_*This*work was prepared under a National Bureau of Standards contract with 
The American University, with the sponsorship of the Army Map Service. 
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An approach to the analytical integration of that equa- 
A new method of fitting the given background curves by 


The storage-outflow discharge relation 
\ general differential equation for water 


0 is derived, with ¢ and 4 constants for the given reservoir, outflow 


The analytical solutions of this 


f(t) are given for the integrable cases (tables 1 
The application of the results obtained is discussed at the end of the 


The applications of eqs (2) and (3) must satisfy 
the basic assumptions made by De Saint Venant in 
their derivation, and specifically that the unsteady 
flow is gradually varied. 

The eq (3) is omitted in many engineering prob- 
lems, as in the routing of floods through great reser- 
voirs and along channels of small slopes with un- 
steady flow very gradually varied, so that only eq (2) 
is used. 

Two important problems which have not as yet 
been completely discussed, and are not considered 
here, are: (a) Under what conditions can eqs (2) 
and (3) be used, and (b) when can eq (3) be neg- 
lected without sensible error. 


The first term of eq (2), 0(Avr)) (Or) = (OQ) /(Oz), 
can, for suitable small Av, be replaced by (AQ) /(Az); 
for a given time interval Af and a reach of length Az, 
can be replaced by AQ/Ar=(P—)/Ar, where P is 
arrival or inflow discharge and @ is the departure or 
outflow discharge in m*/s or in cfs. (Here P and 
() are used instead of usual J and 0.) The second 
term (0.4)/(Of) can be replaced for small At by (4A)/ 


(At)=(AA,+ SA,)/(2At), where AA,=A;’—A; and 
AA,= A}’— Aj, and AA, and AA, are changes of 


cross-sectional areas at the beginning, and at the end 
of the reach Ary during time Aft, and A; and A) are 
areas for the beginning, and A,’ and A,’ areas for the 
end of time interval Af. Equation (2) can thus be 
expressed as 


p—qai4 





tA. 1, Ar_AS 
2 At At’ 


where AS is the volume change for the reach Az 
during time interval Af. If now Aft--0, this differ- 








ential equation becomes eq (1), in which S is the 
storage volume for body of water of the considered 


reach Az, and dS/dt is the rate of change with time of | 


this volume. In the case of a reservoir or lake, with 
small velocities along the lake and with inflow com- 
ing from different tributaries, the reach Av is replaced 
by the water body of the lake. Strictly speaking, 
some end portions of lakes with shallow water and 
high inflow should be considered separately, but their 
influence on the outflow hydrograph can be neglected 
in comparison with the main influence of great 
fluctuating storage space. Equations (2) and (1) 
are identical, and the equation of continuity in its 
form of eq (1) will be called the sterage differential 
equation for water in lakes and along channels of 
small slopes with unsteady flow very gradually 
varied with time. 

The analytical integration of the storage differ- 
ential equation would be useful in the solution of 
problems in cases where its terms can be approxi- 
mated by integrable expressions. Some of these 
cases are: (a) Routing of very gradually varied floods 
along large channels of small bottom slope; (b) effect 
of unregulated lakes on flood waves; (c) routing of 
flood peaks through reservoirs for the inflow dis- 
charges higher than that discharge of hydrograph, 
which corresponds to the full capacity of open gates 
and valves; (d) study of the genesis of hydrograph 
and separation of water flow according to its origin 
(surface runoff, underground storage, lakes, channel 
storage, ete.); (e) outflow through partial openings 
in case of rupture of dam with no influence of tail- 


race levels on the outflow hydrograph, when the wave | 


movement along the lake created by rupture can be 
neglected; (f) computation of the seepage water out 
of reservoirs, ete. 

Equation (1) serves generally for the computation 
of relations between five functions: 

(1) Inflow hydrograph, P=/f, (1); 

(2) outflow hydrograph, Q=/f.(?) ; 

(3) stage hydrograph, //=f,(¢) ; 

(4) outflow rating curve, Q=f,(/7); 

(5) storage function, S=/;(/7); or area function 
A=f,(11), with five variables, Q, P, /7, S, t. 

When three of five functions with boundary condi- 
tions are known (three variables can be excluded), 
the eq (1) enables the computation of the relation 
between two remaining variables. 

The usual integration procedures are numerical, 
graphical, combined numerical and graphical, and 
integration by special devices. 

The analytical integration of eq (1) is net usually 
feasible in practice, mostly because of the difficulties 
of fitting easily three known curves out of five by 
tractable and integrable mathematical expressions. 

The most common cases in the application of eq 
(1) have, as given, the following curves: (1) Storage 
or area function (obtained by survey); (2) outflow 
rating curve (obtained by gauging, hydraulic com- 
putation, model study, ete.); and (3) inflow (or out- 
flow) hvdrograph. Two other functions are normally 
to be computed; (1) outflow (or inflow) hydrograph, 
and (2) stage hvdrograph. 


44 


The subject of this paper is the analytical integra. 
tion of eq (1) for those cases for which inflow hydro. 
graph, storage function and outflow rating curve egy 
be entirely or partly fitted by simple expressions 
which make eq (1) integrable. 


2. Fitting of Mathematical Expressions to 
Given Curves 


2.1. Storage Function 


The storage function, which relates lake or channe! 
volume to its level referred to some datum, can be 
approximated either by the function (see references 
[1] ' to [5]) of the type 


S=all"™ (4) 
or by the polynomial of the type 
S=Ap+A,H+ AIP 4 +- A,,f1™ (4a) 


where S=storage volume, //=depth of water above 
a reference level suitably selected, a@ and m, and 
Apo, A, As, oe A. are coefficients to be determined 
from data of storage function. 

The reference level of eq (4) is normally that of 
zero storage (lowest level of lake or reservoir, or 
river bed, ete.). The level of zero outflow does not 
coincide generally with the level of zero storage, and 
in that case another form of eq (4) is used: 


S=a(H"—Hy?) (4b) 
the level of zero outflow 
and of the level of zero storage. The coefficients a 
and m in eq (4b) are the same as in eq (4). The 
use of eq (4), as will be shown later, gives a form to 
eq (1) which is integrable in many cases, but eq 
(4b) makes it less simple and more complex for 
analytical treatment. 

Equation (4a) is well-suited to be used as the 
storage function, if enough terms (m-+1) are in- 


where //,—difference of 


volved. It must be supposed that for the level 
/1/=0, the storage is 4p, so that eq (4a) can be used 


in case the level of zero outflow coincides with the 
level of storage Ay. If the level of zero outflow is 
changed, the coefficients Ay, A;, Ae, .. A, have 
also to be changed. 

Instead of using eq (4b) repeatedly with various 
values of //, to cover the entire range of S versus H, 
a family of eqs (4) may be used with continuously 
changing a and m as function either of level of zero 
outflow or of //,, the distance from level of zero 
outflow to a reference level. The coefficients a and 
m of eq (4) depend on level of zero outflow, or on 
difference of levels /7,, and are determined by the 
condition that /7=0 of new zero storage coincides 
always with the level of zero outflow. 
that the coefficients a and m are to be determined 
for different levels of zero outflow (and therefore also 
of zero storage) as continuous funetions. The 


Figures in brackets indicate the literature references at the end of this paper. 
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question arises Whether eq (4) fits the thus conceived 
changing storage function as well as it fits the entire 
function from the lowest to the highest level of the 
body of water. In some cases eq (4), used in the 


manner described above, fits the upper parts of the 
storage curve better than the entire curve, but in 
other cases the reverse is true. 

Figures | and 2 show the second case, and figures 
9 and 4 the first case. 
functions S 


Figures 1 and 3 show the 
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Figure 2. Storage function of figure 1 in logarithmic scales for 
different values H, of zero outflow, where the storage functions 
from H 
S=aH» (first excam ple ) 


up to the highest level are fitted by the power functions | 
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functions m 


J (H.) and a 


depth of water in reservoir at dam, and the 


f.(H.) with H, as dis- 


tance of level of outflow zero to the sea level, or to 


the reservoir bottom at dam /7{. Figures 2 


and 4 


show the fitting of eq (4) to the parts of storage 


curves for different values /7/,. 


When //, is ex- 


pressed as //;, or as distance to the bottom of reser- 
voi, it can be used as dimensionless number h= 


H/H,,, with H,, 


normal operation. 
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Storage function of figure 3 in logaiithmic scales for 
different values H, of level of zero outflow, where the storage 
functions from H, up te the highest level are fitted by the power 








The storage functions in figures (1) and (3) are 
given from the reservoir bottom at dam up, and in 
figures 2 and 4 from the level of zero outflow up 
for each of given curves (as straight lines). 

The funetion m=/,(/7.) decreases with higher 
value of /7,, i.e., higher level of zero outflow, and is 
approximately 1 for the highest levels and small 
range of level differences, where the storage function 
is nearly a straight line. For lower //, and for higher 
range of storage levels, the value m is greater. The 
function a=/f,(/7.) has the opposite characteristics. 

If only one value of //, is to be investigated, a and 
m are fixed coefficients, but if a range of //, is to be 
investigated (dam-breaching problems with various 
breachings and different lowest outflow levels), the 
functions m and a of 7/7, can be easily obtained as in 
the case of figures 1 to 4. 

When the sedimentation of reservoirs has to be 
taken into account, a and m are not only a function of 
level of zero outflow, but also of time (apart from 
sediment characteristics they depend on operation 
practice of reservoir, inflow hydrographs, etc.). As 
a and m make eq (4) less tractable in case they change 
with time (which is sometimes rather difficult to 
predict), they will be considered here as constants for 
a given reservoir, for given level of zero outflow, and 
for selected time period of reservoir life. 

The coefficient m ranges within the limits 1 to 5 
but for the majority of natural valleys the range is 
about 1 to 4. When //, is high and the range of 
levels is very small, m is mostly 1.0 to 1.5 and rarely 
greater than 2; for highest range it is 2 to 5. The m 
depends not only on the range of levels, but also on 
the shape of reservoir cross sections. When those 
sections can be fitted by a power function for width 
of type B=2pIl, m is function of s. The value a 
depends on //,, but apart from that, @ is higher for 
wider valleys, for lower river slopes, and for lower 
values of exponent s. 

The property of the storage function to be well- 
fitted by an expression of the type of eq (4) can be 
used to compute this function in the case where only 
two points of the relation storage versus level are 
available. These two points can be found in the 
literature for practically any reservoir: capacities at 
highest and at lowest operational levels. Two pairs 
of values (S, H/) are sufficient for the computation 
of values a and m in eq (4), thus allowing the deter- 
mination of the curves of figures 1 and 3. Taking 
into account the errors in the survey of reservoirs 
and the fact that the storage function is being con- 
stantly changed by sedimentation and littoral erosion, 
the values @ and m determined by eq (4) and the 
figures 1 to 4 are accurate enough for the analytical 
solution of practical problems. 


2.2. Outflow Rating Curve 


The outflow rating curves are relations of the 
departure discharge to the water level in the res- 
ervoir. They can be fitted also in many cases 
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either by an expression of type 


Q—b/T’, 5 
or by a polynomial, 


= by +b 1+ bSP +b fT -h HH’, da 


where Q=departure discharge, /7—level above thp 
depth for which Q=b), b and r or bo, b:, bo, bs, 
b, are coefficients depending only on boundary 
= (outlet shape) and on type of flow (closed 
‘free outlets). The eq (Sa) can be also conceived 
as a rating function expressed in power series form. 
The range of r in eq (5) for free surface outflow js 
1.5—4.0 (usually 1.5 to 3.0) and for pressure type 


flow about r=0.5 in most cases. The range r= 
0.5—1.5 can be covered also in spec ‘ial cases for 


special shapes of outflows. 

In case two types of flow, 
and free surface flow, occur 
takes the form 


viz, flow under pressure 
simultaneously, eq (5) 


6,71," + 6.97," =6, (H+-H,)"+6b.H", (db 


4 
where /1,=H and /7, 
heads (level differences) to be 
and latter types of flow respectively, 
correspond to closed pressure flow and 6, 
free surface flow. As 7, and 7, are not 
integers, eq (Sb) is different from eq (Sa). 

The use of eqs (5a) and (5b) instead 

makes eq (1) integrable in a much smaller number 
of Equation (5) will be used exclusively 
here. 


I1+-H/, are the hydraulic 
taken for the former 
b, and +, 
and ry to 
generally 
of eq (5 


Cases. 


2.3. Inflow Hydrograph 


It is practically impossible to fit an entire natural 


inflow hydrograph by a single mathematical expres- 
sion, but it is possible to fit some of their parts by 
tractable and simple functions. 

The simplest forms of function P=/(t) to be used 
for integration are: ) PO (arrival discharge zero 
as dry period of year or as water fully stored in } 
upstream reservoirs, etc.); (b) P=Py=constant (or 
nearly constant, with low river flow slowly changing, 
regulated constant flow from upstream reservoirs, 
short-term operation of reservoirs, when flow could 
be placed by constant discharge); (¢) P= P)—ft 
(gradually varied flow); P= Po Pura": # 
Pit'e-" (Po, f, and s are different constants in each 
case ). 


2.4. Storage-Outflow Discharge Function 


When the lowest outflow level coincides with the 
level of zero storage, or when reference levels fot 
measuring /7 in both eqs (4) and (5) are the same 
(adjustment to this condition is always — by 
taking the corresponding coefficients a and m from 
their functions, fig. 1 to 4), the dimimation of H 
gives for the storage-outflow discharge function, or 
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where c= 6°" /a and n=m/r. 

The coefficient » for free surface flow is in the 
range 4 to 4, and for closed pressure flow 2 to 10, 
or a little greater. 

The use of pairs of eqs (4) and (Sa), (4) and (5b), 
ja) and (5a), (4a) and (5b), (4b) and (5), (4b) and 
a), (4b) and (5b) gives a more complex relation 
s=/(Q), in most cases in parametric form with // 
as the parameter, than the use of (4) and (5) or (4a) 
and (5). However, in this last case 


| I 
S ay+- ()' rt 
C; C» 


(Oa) 


l 
.+ 


Cc 


"+... | al 


m 


which makes eq (1) less tractable than the use of 
eqs (4) and (5). 

As eqs (4) and (5) fit well many storage and out- 
flow functions, for different value of /7., and 
eq (6) allows the integration of eq (1) in many cases, 
eqs (4) and (5), and (6) will be used in this study. 


as 


3. General Type of Differential Equations 


With the introduction of eq (6) into eq (1), one 


obtains n Q"~' dQ+e(Q—P) dt=0. 
The substitution y=Q-" and the replacement 
k=(2n—1)/n va. =2—— then gives the following 
n mM 
expression: 
y’ +ePy?—cy* —0. (7) 


This is the general differential equation for storage 
reservoirs and storage channels with no artificial 
control of outlet flow, under conditions discussed 
previously in this paper. For free surface outflows 
the range of & is usually from about —2 for n= to 
k=", for n—4, and for the closed type outflows the 
range of k is from *% to ' > or a little greater, but 
less than the limit value k=2. Practically, the range 
of k is from —2 to +2, where k= -+-2 is not possible. 

The relations of k or » to m and r are given in 


figure 5. The most common range for both m 
| to 4) and r (0.5 and 1.5 to 3) is specially empha- 
sized for reservoirs in river valleys. In most 
practical cases k>O, and usually k>>1/2. 
The eq (7) could also be expressed as 
y’ +ey? (P—1)—ey'/"=0. (S) 
The general form of eqs (7) and (8) is 
y’ + AM + NyF¥=0 (9) 
Where J and N are functions of ¢; 4 and g are any 
real numbers. By the substitution y=Z7°/">" the 
general form is reduced to 
h+g—2 
7! + (h—-1)MZ?+-(h—1)NZ *-' =0 (10) 
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Figure 5. Values n of storage-outflow discharge relationship 


, in elation to m and r (exrpo- 
n 


9 


i 
s Q", and values k 
c 


nents of storage and of rating curve, expressed in power form). 


which is eq (7) under conditions that A—1 
N=-—1, M=P, and (h+g—2)/(h—1) =k. 

The general analytical solution of the ordinary 
differential equations of the type (7) or (10) for any 
k and for any expression of P does not exist. 
Chiellini [6] has shown that the solution is possible 
for k=3 if the expression (\M’N—AMN’) is different 
from \/° by a constant factor. When & is a positive 
and integer number higher than 3, Luigi Conte 
(6, 7] has given the solution for special cases. As & 
is always less than two here, their’ procedures can- 
not be used for eq (7). 

Both facets, (a) that & is a number derived from 
continuously changing natural conditions and is 
rarely an integer, and (b) that the function P=/f(t) 
is not simple and must be fitted by integrable 
expressions, make the analytical integration of eq 
(7) possible only in special cases. For k=—2, —1, 
0, and 1 and some rational numbers the integration 
is simpler for some types of P, but the integration 
for intermediate values of k and for the majority 
of P functions has to be performed numerically or 
graphically. The analytical integration for special 
values of & and for selected functions P=f(t) is 
useful, because the resulting functions give the type 
of curves which can be expected for intermediate 
values of k&. Or, when & and outflow functions are 
known, some conclusions can be made about the 
inflow functions. The type of functions obtained 
in integrable cases can serve for control of curves 
obtained by other procedures. 


=¢C, 





4. Integration in Special Cases Q/Qy depends on both ¢ and Qo, except for n=k~ 


(table 1), in which case it depends only on ¢. 10 
Equation (7) will be analyzed for some values of The analytical expressions of eq (12) for some } 
k and some types of the function P?. The values | (Or ”) values are given in table 1 . Figure 6 gives 
of & taken will be 2, —1, 0, 1 and, tentatively, | curves for n=}, 1, 3, 2 and 3 for given values of os 
—3/2, —1/2, 1/2, 3/2 along with some other rational | @ and 6, with e=6"/a. The y ‘alues a of eq (4) and} 
numbers. | of « eq (5) are taken constant for all five curves (fo 
41. P—0 easier computation in examples), though they change - 


when m and r (and so n) changes. For n 2 the 
outflow hydrograph is a straight line, which corre. 
sponds, among other cases, to the cylindric reservoir 
, kg. (11) | (m=1) and to the outflow under pressure (r=}) 04 


Equation (7) becomes 


y’ —cy 
The same type of linear hydrograph is obtained for 
with y=y=Q”" for t=0, where (, is outflow dis- | “%=3 and r= 3, which is the case for some deep 
charge at the time t=0. The solution of eq (11) | lakes with outflow function Q=6b/72. For river oz 
with y=Q-" is then outlets from lakes, with r>$ and m<3, n<2, lj 
hydrographs are concave upwards. For n=1, the 
Qr'=@r i ;—! ) et (12) solution of eq (11) is an exponential function. mm) 
0 n : " functions for n <1 are asymptotic with the time axis 
Hvdrographs for »>1 do not have this asymptote, — 4.: 
This is the general expression for outflow hydro- | but cut ker" —axis for finite | (for n=2,5 | 
graphs from storage reservoir with no inflow into the | f=2@o/¢). _ For the range 2 >n the hydrographs»: 


reservoir. The solution can be obtained from eq | are power fune oon and concave panes but have 
(12) for each value of 7 except for n l (k—=1). finite values for QO (are not asvinptotic to the? 
but eq (11) gives solution in this case directly | ®XIs). For n<l an curves are hyperbolic functions, — ,. 
(table 1). Therefore, for P=0, or when the constant | asymptotic and concave upward. 

All outflow hydrographs for very small water levels °° 
above the bottom of outflow orifice become asymp- Fic 
a—!l a y 9 ial totic, because their outflow functions change due to ( 


— * % (1 n (Qi the change of type of flow or to some secondary } ou 
effects. In closed- type outflow the function changes 
. “ ° - | from r=} to about r=% as soon as the water | wh 
level drops to near the upper edge of the outflow eqs 
2 4 ol 14800 41) ~3 orifice, and changes further for smaller outflow | as: 
, : alssect)-3 heights due to surface tension and other secondary | inv 
=e effects. The exponent r increases so that n= m/r<1. l 
; 0( 144 0 ‘r) ~§ The shapes of outflow hydrographs depend thus . has 
5 on the ratio of exponents of storage and outflow nu 
, oe Seen © funetions, and they ean be: hyperbolic (n<1), 
; exponential (n=1), power function with exponents 
1 ; do( 14 S0.b¢) higher than 1 (2>n>1), straight line (m2) and 
‘ power function with exponents lower than 1 (n>2). I 
, Owe ~t They show that the natural conditions can produce int 
bo tL different outflow hydrographs on the same river, 
o{ 1-—! ) depending on the lake and channel characteristics 


and outflow shapes. 

The free spillways of reservoirs, or free outflows — anc 
of shallow lakes, where the range of storage fluetu- 
ation is small, m=1.0 to 1.3, and r> 3, have n<l, 


nh 
~ 
— —— — a. 
te 
~ 
— 


2c 5 : 
' a 1-375!) and all outflow hydrographs are hyperbolic asymp- 
totic functions (fig. 6). For the higher range of 
ae 1 z : : ? ; 
7 ' of 1 ‘) levels of outflow, m is greater, and the outflow ” wit 
10, »* : 
hvdrographs are power functions whose graphs are y 
; ; . convex upward or downward. The smallest value of | car 
inflow discharge is small and can be neglected in ' | the 
: : : = ; mee 9 ‘ ‘ 
comparison to the outflow discharges, the analytical |” for closed-type outflow is x Sith baad that af ana 
ri ako of rk, 4g ow solutions for —_ k(n), hvdrogri aphs are powe r functions convex upwe ard, IS U 
“oo o "ion | P ‘ l f In case eq (Ga) is used instead of eq (6), the follow- ; rau 
quation (12) can De written in the form: ing general solution of eq (7) for P=0 is obtained: Inte 
| i ft 
nN 1 eft eal ‘ 0 
Q Qo ( — (i; } (13) m— . 
" 0 ~~ m r ’ 
f-= > a Q , (14 


This shows that the decrease of discharge ratio 
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FigtrE 6. Outflow hydrographs in their recession parts from 


(Q), to zero, related to different values k and n, in the case where 
inflow is zero 


which for the case m> 3 is more complicated than 
eqs (12) or (13), particularly since eq (14) gives ¢ 
as a function of Q@ rather than the more convenient 
inverse. 

If both types of flow occur simultaneously, eq (5b) 
has analytical solutions in a much more limited 
number of cases, 


4.2. P=P,—Constant 


In the case of constant or nearly constant inflow 
into the reservoir, with P=P,, eq (7) becomes 


y’ +ePyy—cey =9, (15) 
and the separation of variables gives 
" @e . 
5 at ptet=K (16) 


with A=integration constant. 

The indicated integration of eq (16) cannot be 
carried out analytically for arbitrary &, contrary to 
the case of eq (11). Equation (16) has, however, 
analytical solutions for many rational numbers, which 
is useful, since it is possible to approximate & by a 
rational number of type g/h, where g and h are 
integers. Equation (16) can be written in the form 
of the binomial integral 


y~*(Po K. (17) 


y~*) dy +et 
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The binomial y~*(P)—y*~?)*! is of the form 


y*(ay’ +b)? (18) 


which can be reduced to an algebraic function in case 
any one of the three numbers 


s+] s+] 
Pp, ——» and +p 
= 
is an integer. As p=—1 is an integer, when & is a 


rational number, eq (16) can always be reduced to 
an algebraic function. 

The most general solutions of the binomial integral, 
if the obtained algebraic functions are integrable, 
are in the form of sum of logarithmic, algebraic and 
arctan terms. In some cases the arctan term does 
not appear in the solution, as is the case also with the 
algebraic term, but the logarithmic term is always 
present. The solutions of eq (16) will be analyzed 
here for some values of k. 

By substitution of a new variable uw in eq (16), 


1 
u=y'*.P,*-9, (19) 
the following general binomial expression of eq (16) 
for k=g/h is obtained: 


g—h 


MP | u®-*&(y*-8§—1)—"'du+et=K (20) 


with h—g—1 and 2h—g integers. 


For 

f(u) | uh #—1(y?*—8#_1)— "dy (21) 

the constant AK for t=0 (u=up) isK=hP,S-”(%-® X 
f(u), and the general solution has the form 

t= Rf (m)—f(u)] (22) 


h ; 
with R=-P,*-”/*-®, 


The following values and ranges of g and h are 
analyzed: (1) g=2h; (2) AS g<2h; (3) g<h; (4) 
g=0, h¥0. 

(1) g=2h. For k=2, eq (20) has no practical 
significance. 


(2) h<qg<2h. The value s=2h—g is always a 
pesitive integer number, and j7=Ah—g—1 is always a 
negative integer number or zero. Equation (20) is 
integrable if s=2h—g<4 or s=6, when the rational 
function can be separated in integrable terms, what- 
ever the value of ;=h—g—1. The solution is 


h bs 
t—— Pp" 


c 


Lf Uo) (22a) 


f(u)). 








(a)s=2h—g=1, n=h. 
, u—l1,1 l l 
— ‘ W1s 
f(uj=In ; -$—+-——+...+ ——- (21a) 
l uo 2u (n—1)u 
with u= P,Q". 
This case covers the following values of » and k: 
Seis & & & i) 7, ete. 
k: 1, 3/2, 5/3, 7/4, 9/5, 11/6, 13/7, ete. 
Some of these cases are given in table 2. 
For n=1 (k=1) only the logarithmic part of eq 
(21a) is present and eq (22a) gives: 
(J 0 T (to Poe = (21b) 
(b) s=2h—g=2, n=h/2. 
u—1,2, 2 2 
f(u)=In - -+——.+ ... —__. 
’ us ity at T3in LT eit 
(21¢) 


with u=P)?Q"'". This case covers the values: 


n: 3/2, 5/2, 7/2, 9/2, ete. 
k: 4/3, 8/5, 12/7, 16/9, ete. 


Some of these are given in table 2. 


TaBLE 2. P= P,= Constant t= R {flu f(n)] 
7) n P 
1 1 1 l 
0. 25 2 Ps P*Q 4 In +2 uretan 
4 ‘ l 
l l 1 l l 
0.33 | P 3 P3¢ 3 in f arctan 
3 ri + : ( ) 
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2 l 23 } F . , l 2u+l 
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(c) s=2h—g=3, n=h/3. The solutions for thre. 
values of h and n (h=4, n=4/3; h=5, n=5» 
h=7, n=7/3) are given in table 2. _ 

(d) s=2h—g=4, n=h/4. For h=5, n=5/4 th, 
solution is given in table 2. 

(3) g<h. The analytical imtegration ean}, 
earried out when 4<2h—g=s<4, or g/2—2<)) ™ 
<g/2+2. For h>g/2+2 and h<g/2—2 the ex. sol 
ponent s=2h—g in eq (20) is higher than 4, so tha 
the algebraic function u®-"*'e"-*—1) can not by 
separated into integrable components except. jj 
special eases. The solutions for n=—3/4 (& 9/2 
TD) 2/3 (k ] Zi. n 1/3 (k 1) and 7] l4 
(hk 2) are given in table 2. re 

(4) g=0, h=1; n=1/2 (k=0). = The solution js 
given in table 2. Figure 7 gives five curves for thy 
values n: 1/2, 1,3/2, 2 and 3, as in figure 6. — All fiy, 
are asymptotic to the value ?)—0.20, but for t neare 
to zero they approach the curves of figure 6. 

4.3. P=f(t) 

The solution of eq (3) for four tvpes of P=fi 
can be obtained for special values of k, but th 
general analytical solution such as was obtained j 
the case of eq (13), or in many special cases of e 
(22), is not possible. 

P= P=constant ; 
Q Q= 1.0 D*0.05 metric 
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The 


a. k=1 (n=)) 
Equation (7 becomes 
y’ +ePy—cy=0, (23) 


vhich is a Bernoulli equation, and has the general | 


solution for y= l. 
( (K Cc | Pr ‘dt ). (24) 
The solutions for four types of function P=f(t) 
are cviven in table 3. 
TABLE 3. P=f 
} ] n l « integ: 
/ fit , } 
p= pP ( F ) +( l , ). 
i } [ « / 1 “] ~ [: _ ] 
f cl 
( yentt ( —) 
[ “re 
Suc! analytical solutions can be obtained under 


condition that the integration of 


| Peé'dt 25 
eq (23) can be earried out in closed form. 
bk 0 (n-5) 
2 
Equation (7) becomes 
y’ +cPy?—e=0 (26) 


Which is a special type of the Riceati equation, with 
y=", As two particular integrals of eq (26) are 
KNOWN, VIZ, Yom Qs‘? for t=0, and Q=0 or Q=P, 
lor ¢ infinite, the solution of eq (26) is possible. 
The substitution 7 Y¥— Yo gives 

ZL! + PRP yZ- 


yo ePyi—e) =0 


and eq (26) becomes the Bernoulli equation 


Z'-+-cPZ2*+ 


i 


2¢PyoZ=0 


with the solution 


dt) (28) 


constant AC is to be determined from the known 
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value of Q for ¢ infinite. 


The solution for P=Py,e~" is 
ji : 
es L\ 
A Or . . (29) 
do oe ft \ | 


where 


q 2el o fy (Wo. 


The analytical integration is not possible for the 
three other P functions of table 3, because the integral 
Se dt can be obtained, in general, only in power 
series form. 

This analysis shows that the analytical integration 
for simple expressions P=f(t) can be done for k=1 
and for /=0 in some cases. 


5. Discussion 


The analysis of the feasibility of fitting the given 
curves, as parts of differential equation for water 
storage, by mathematical expressions and its inte- 
gration by analytical procedures in special cases, 
show some characteristics of outflow hydrographs. 

The shape of outflow hydrographs from lakes and 
channel reaches is highly influenced by the ratio of 
exponents of storage function and outflow rating 
curve, both of the power type y= sz’. 

The analytical integration of the storage differen- 
tial equation can be performed for any ratio of ex- 
ponents in case of zero inflow, and for many ratios 
of exponents in case of constant inflow. The inte- 
gration is, however, limited to a small number of 
ratios and inflow functions, when the inflow changes 
with time. 

The usual procedure, to plot the recession part of 
river-flow hydrographs on semi-log paper and to fit 
it by a straight line, shows that in the majority of 
cases the fitting of a straight line is a rough approxi- 
mation. The results of above analytical mtegration 
show that the straight line can be used for accurate 
fitting only if the ratio of exponents of storage func- 
tion and outflow rating curve has a mean value near 
unity (n=m/r=1) and for zero or constant inflow. 
This can be considered a special case. 

Though the discharge hydrographs of underground 
water in a river basin can be (theoretically) fitted in 
many cases during the recession period of no water 
supply to the underground by an exponential fune- 
tion of type Q,e7"', the water storage either in lakes 
or along river channels influences highly the shape 
of river hydrographs, so that it becomes less and less 
of the pure exponential type. The higher the ratio 
between the effective storage (storage which influ- 
ences the outflow hydrograph) of lakes and channels 
to the total river flow during the period of recession 
flow, and the higher the departure of the ratio of 
two exponents from unity, the greater is their influ- 
ence on the shape of river hydrographs and the larger 
is the departure from the original hydrograph of the 
recession flow from the underground. 

As the ratio n=m/r changes along the river 
channel, the accurate flood routing procedures must 








take that fact into consideration. 


The assumptions usually taken for some breach 
. S| 


The described analytical procedure is useful for the | and for inflow discharge, and the accuracy of basi 


study of outflow hydrographs in case of rapid 
openings in water bodies (breaches of dams, rapid 
openings of gates and valves, breaches of channel 
walls or levees, etc.). The assumption of simple 
shape of openings and of simple inflow hydrographs 
(zero, constant inflow) usual in case of dam 
breaches. As the study of dam breaches has to be 
made for many openings of quite different dimensions 
and shapes for different inflow discharges, the 
analytical integration as given above can have some 
advantages in the computation of hydrographs in 
comparison to the standard graphical or numerical 
procedures of integration. 

The difficulties in fitting the background curves 
by tractable mathematical expressions and the 
difficulties of integrating analytically the storage 
differential equation, and especially in ease the inflow 
changes with time, explain why the usual graphical 
and numerical procedures have taken so common a 
place in the engineering practice. Nevertheless, the 
analytical treatment of the storage differential 
equation can be useful for some problems, can serve 
for better understanding of relations, and can save 
efforts in problems, as, for instance, in case of outflow 
hydrographs out of dam breaches. 


Is 
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data in the routing of waves createa downstrean 
justify the use of the procedure given for the fitting 
of mathematical expressions to the basic and give; 
curves and the use also of the analytical integratigy 
of the storage differential equation for the com. 
putation of outflow hydrographs. 
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Compressible Turbulent Boundary Layers With 
Heat Transfer and Pressure Gradient in Flow 
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The best-known theoretical works on boundary-layer problems, especially in the case 
of compressible flow without or with heat transfer, are related to the laminar boundary 
laver, although the turbulent boundary layer is, in practice, often more interesting. The 
laminar boundary layer is more easily accessible to theoretical treatment because clearly 
defined relations exist between the viscosity w and the shear stresses r. In the turbulent 
case, empirical relations must be introduced. Therefore, attempts to get exact solutions 
are not worthwhile, while efforts to obtain approximate solutions, based for instance on 
the von Karman-Pohlhausen principle (Z. Angew. Math. u. Mech. 1, 233 and 252, 1921) 
of utilizing integral conditions, appear to be appropriate to this problem. In the last few 
years the accuracy of such approximate solutions for the incompressible case was notice- 
ably improved by the application of a new energy integral condition in connection with 
a new empirical law for the dissipation in turbulent boundary layers, stated by J. Rotta 
(Ingr-Arch. 20, 195, 1952) and E. Truckenbrodt (Ingr-Arch. 20, 212, 1952). The empirical 
laws for dissipation and for turbulent wall friction, which are needed in the present approxi- 
mate theory, are formulated on the basis of available measurements for incompressible flow. 
Generalization to the compressible flow with heat transfer is made from physical considera- 
tions. Calculated results agree satisfactorily with available experimental data. Some pos- 
sibilities for improving as well as for simplifying the approximation theory are outlined. 


1. Introduction 


In principle, problems of compressible /aminar boundary layers may be solved by using 
the partial differential equations for the velocity and temperature fields, and the thermodynamic 
relations between viscosity u, density p, pressure p, and temperature 7. Exact solutions for 
special cases have been worked out by several authors [5, 6, 7, 8].2 Nevertheless, approximate 
theories of the von Karman-Pohlhausen type have been successful because of their easy adapt- 
ability to the most general case of flows with arbitrary pressure gradient in flow direction 
{9 through 16). 

An example is the case of a compressible turbulent boundary layer, which is more important 
in practice because the controlling viscosity is no longer a material property, such as the molec- 
ular viscosity, w; rather, it is like so-called “apparent viscosity,” g,, which is essentially 
determined by the empirical laws of turbulent mixing motions. In view of the empirical 
character of these laws, attempts of exact solutions of Prandtl’s boundary layer equations 
are indeed inappropriate. Hence, approximation theories of the von Karman-Pohlhausen 
type are indeed the only ones which are well adapted to the problem of calculating turbulent 
boundary lavers. 

The first solution was given by E. Gruschwitz [17] in 1931. He introduced a one-parameter 
family of empirical turbulent velocity profiles with parameter /7, and, as the second param- 
eter of the problem, the boundary-layer momentum-loss thickness 6. Gruschwitz then 
derived an empirical energy law. This was used together with the known von Karman integral 
condition for the momentum within the boundary layer, to determine the two parameters 
IT and 6 as functions of x for a given pressure gradient. 

The so-called Gruschwitz energy law could not be verified in a general manner and appears 
therefore somewhat doubtful. Moreover, it is restricted to the incompressible case. A definite 
improvement of the approximate theory for turbulent boundary layers was accomplished 


This work was supported by the U.S. Air Force, through the Office of Scientific Research of the Air Research and Development Command, 
2? Figures in brackets indicate the literature references at the end of this paper 
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when several authors [9 through 15] independently (between 1935 and 1946) showed that 
an integral condition of the energy within the boundary laver may be derived from the Prandtl 
boundary-layer equation in the similar way as for the momentum. Some years later (1950) 
J. Rotta [3] derived from measurements the empirical law for the dissipation within the turbulent 
boundary laver, which was needed for the applications of the new integral energy condition. 
At about the same time the empirical law for the turbulent wall friction occurring in the mo- 
mentum condition was improved by new accurate measurements of Ludwieg and Tillmann 
[18], such that it may be generally valid for positive and negative pressure gradients in flow 
direction. 

Now, on the basis of the momentum and energy integral conditions on the one hand 
and the two empirical laws for the turbulent wall friction and dissipation on the other hand, 
a two-variable method approximate theory of boundary laver could be developed, similar 
to that formerly worked out by several authors for laminar boundary lavers |9 through 16). 

For the incompressible case, the integration method of E. Truckenbrodt [4] (1952) (similar 
to that developed by A. Walz [19] in 1941 for the laminar case) appears to be the simplest. 
However, applications to compressible turbulent boundary laver failed, because the semi- 
empirical laws for wall friction and dissipation could not be set up, since before that time 
measurements of turbulent boundary lavers at high Mach numbers were not available in 
sufficient generality, 

Basing his analysis upon simple physical facts about the turbulence mechanism the 
author was able to derive a generalization of the empirical laws for wall friction and dissipation 
to include compressibility, [26, 27]. 

The inclusion of the case of heat transfer, however, implies a new difficulty msofar as the 
integral energy condition derived from Prandtl’s equation considers only an energy flux along 
the body surface but does not relate to an energy flux through the wall. In the latter general 
case we must start from the general differential equation for the temperature field. In addition 
it must be checked whether the above compressible considerations would include heat transfer 
at the wall. 

As suggested by L. E. Kalikhman [13] for the case of laminar boundary layers the partial 
differential equation for the temperature field can be transformed by partial integrations into 
an infinite svstem of ordinary differential equations similar to that obtained by partial integra- 
tion of Prandtl’s boundary laver equation. 

Hlenee, two infin ite systems of ordinary di there ntial CQ Uations with a twice infinite number of 
unknowns are available which appear to he complete ly equivale nt to the two partial differential 
equations for the velocity and temperature field. 

Some general properties of these svstems, which had not vet been outlined in the literature, 
will be discussed. 

As a first step of evaluation of the possibilities involved in these systems of equations, a 
“two-equation, two-unknown” method is investigated. This is about as simple as the known 
methods for the incompressible case [4, 19]. Nevertheless results obtained are of sufficient 
accuracy for engineering purposes. 

The application of the new method is illustrated in some examples. Directions for improve- 
ment as well as for further simplifications are indicated. 


2. Partial Integration of Prandti’s Boundary-Layer Equation for a Compressible 
Two-Dimensional Steady State Flow Along an Insulated Wall 


2.1. Derivation of an Infinite System of Ordinary Differential Equations 


As proposed by Leibenson [9], Golubew [10], Sutton [11], kK. Weighardt [12], and E. Truck- 
enbrodt [4], Prandtl’s svstem of partial differential equations for a compressible two-dimen- 
sional flow along an insulated wall may be transformed into an infinite system of ordinary 
differential equations by partial integrations. We will give here a new derivation of this 
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system which is appropriate for our purposes. 

Let rs, y be the coordinates, the z-direction being taken everywhere tangent to the wall. 
Further « and ¢ are the components of the velocity vector, respecting in the z and y directions, 
uw the viscosity coefficient, 6 the density coefficient. Then Prandtl’s boundary-layer equations 
may be written as: 


Ou 
, Ou | Ou dus o( "Oy 1) 
’ ’ Us ; ( 
. "as Oy ome dx oy 
O(pu) , Ol pr) 


0 (continuity equation). (2) 
or | dy on 


Here the subseript 6 denotes that the quantity is referred to the edge of the boundary layer. 
We note that this system is valid even for turbulent boundary lavers, if the velocities uv and ¢ 
as well as pare considered to be the mean values over time and if vi is the effective value =e, 
which consists of the molecular viscosity uw, the apparent viscosity ue. Thus, 


Me—=-MT Me - (3) 


For simplicity, however, we prefer to use in both cases, laminar and turbulent, the same symbols 
u,v, p, and w, knowing that in the turbulent case we have p=up. 
First we multiply eq (1) by wv’ and add to it eq (2) multiplied by w*'/(y+ 1). 
Partial integration from y—0 to y= @ then vields after elementary transformations, 
] d . dusdr 


=P) (psu "f+9, +¢e,=0 (p==O.1.7,..«,8@)2 (4) 
pals ~ ad Us 


Herein the following notations are used: 


"“s/u\’ O T 

€,= (v4 | ( ) ( )du, (2) 
a Us OU \ psu; 

u 8 tlhe : 

pu *( ) dy, (t)) 
Pals ps ry : 

F u u \' - 

4, (v4 | E —( ) Jew, (7) 

PLAT) palls Us ; 


m - , . (S) 
OY 

From the structure of the integral terms (5) through (7) we see that for practical purposes it 
will not be necessary to extend the integration process bevond a certain distance y=6, for 
which we have uw=us, p=ps, 7=—0, because for 1 >6 all integrands became practically zero. 
The laver thickness 6 is somewhat vague in view of the asymptotic behavior of u(y) at the 
outer edge. For approximate theories, however, it is very useful and may be defined as wall 
distance y at which w reaches the value vw. with an accuracy better than 99 percent. Hence 
we pul 


0.9901. < Us< U.. ({)) 


2.2. Discussion of the Infinite System of Equations 


Putting y= 0 we obtain the usual momentum integral condition of von Karman-Pohlhausen, 
for v1 we have Weighardt’s energy integral condition. Only these two cases have been 
considered up to the present, but it is clear that all other remaining integral conditions of the 


mav have also values between ind +, even nonintegers as 0.1, 0.2, . ind so on 
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infinite system (4), the so-called higher momentum conditions (yv=2, 3, 4, ...,2°), would be 
important too, since only the complete system is equivalent to Prandtl’s equations. 

Let us first investigate the meaning of the equations of order v, higher than unity. 

First eq (4) is rewritten in a more convenient form. At the outer edge y= 6 of the boundary 
layer, we may assume adiabatic conditions of state. According to the basic laws of gas dy- 
namics we have between the density of ps, the velocity us, and the Mach number A/5-= us/as 
(as=local sound velocity) the following relations (see, for instance, R. Sauer [20]): 

] dé l dus VE 


ve - Mé?. 10) 
pidr ui dr ° | 


Substitution of (10) into (4) gives 
a ae 
fit+ =f, [2 tp M+" te 0 (11) 


where the prime denotes differentiation with respect tos. 

The integral terms e¢,, f,, g, represent the unknowns of the system. They contain the 
velocity and density distributions u(ry)/us(7) and p(r,y)/ps(4), respectively, which are indeed 
the real unknowns of the problem. 

The only way to obtain a solution of the problem was to introduce approximations for 
u/us and p/ps according to von Karman-Pohlhausen’s idea, for instance by means of a one- 
parameter family of experimental curves. In the incompressible case (p/p;—1) the problem 
was indeed solved in this way by taking two equations from system (11) (for »y=0 and »y=1),' 
using as unknowns the boundary laver thickness 6 or one of the functions e,, f,, g,, and, as the 
second unknown, a shape parameter /7 of the velocity profile. In this case the ratios ¢, 4, 
fv/6, gv/6, as well as ¢,/f, and y,/f,, are only functions of the shape parameter /7. 


The compressible case involves new variables p and T which are related by 


7; 
E. (133) 


since Op/Oy=0. Additional equations will be needed to make the problem determinate. 
This is done for the thermal boundary laver. While we will return to this problem in section 
3 of this report, let us first discuss the principal properties of the infinite system (11). 

The order number v may for instance be increased toward ©. In this case we have from 


equations (6) and (7), 


(f(v+2—Mj) +], S.(1—M}), (14) 


and 


because 


u of - u 
( ) 0. ford ft. (16) 
’ > 


Us 


4 The original PohlIhausen method only used the equation »=0 (momentum equation) from this system, As the second equation, Pohlhausen 
introduced the so-called wall condition, 


which follows directly from Prandtl’s eq (1) for the point y=0. It is obvious that this condition, valid for the wall point, would be doubtful if used 


for fixing the shape parameter. Indeed essentia] improvements of results were observed when this wall! condition was neglected and replaced by 
the equation for y=1 (see, for instance, [16]). 
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whatever the functions w/us and p/ps may be in detail. 
The integral expression (¢,) is found to be zero as follows: 


Ya 


u \? 
[ a+» (. ) | >(. 
us »>o 


Since furthermore the shear stress 7 and its derivative 07/Oy are vanishing at y=6 (u=us) for 
arbitrary velocity profile u/u;, we have 


; u ca 
For 0< —<1 it vields 
Us i 


(e,),»5..—0. (17) 
Hence in the limiting case of y=, system (11) reduces to 


fit (1—M?2) =). (18) 


» . x 
Us 


This equation may be called the mass-flow equation and is nothing less than an averaged form 
of the continuity equation. With the aid of (10) it can immediately be integrated to 


f. pss (6—6,) psus—const. (19) 


where 
“6 pu , 
=| (1-4) dy (20) 
J0 pss . 


is known as the displacement thickness. This relation has not previously been used in approxi- 
mate boundary layer theories. It expresses that the mass flow between the stream lines y=6 
and y=-6, must be constant. We have not ascertained whether relation (19) will be very 
useful in practice. But it is clear that if the velocity and density profiles are once determined 
in the inner part of the boundary laver (thereby determining the value 6,), the total boundary 
laver thickness 6 is determined also by relation (19). 

We have found that if the two-equation system »=0, y=1 with two unknowns is used, this 
relation is mostly violated. Nevertheless the essential properties of the boundary layer, for 
example, wall friction and separation behavior, appear to be correctly described. 

From this consideration about the mass-flow equation as limiting case for y=, we may 
learn that the higher the order v, the less details of the velocity profile will be included in the 
corresponding equation, because the term e¢,, which contains the shearing stress (and with 
this the entire boundary laver effect), approaches zero, with v-o. 

It is therefore not surprising that equations corresponding to low order, for example the 
equation for »=0 and v=1, still give a very good description of the properties of the boundary 
laver if a rather physically reasonable family of velocity profiles is introduced into the equations. 
A good way in principle as suggested for the incompressible laminar boundary layers a long time 
ago by several authors [9,10,11,12], would be to restrict to a finite number of equations of the 
system, preferably to those of low order, say up to v=5 or 10, and to introduce a polynomial 
expression for the velocity profile 


M- ral YH (xr) ,Ho(r),H;(2)...., H0) | (21) 


Us Us Lbo(r) 


with » coefficients 7. Then »+1 ordinary differential equations will allow the determination 
of the » coefficients 77, and 6 as functions of x. In this case no detailed assumption about the 
velocity profile, except some trivial boundary conditions, must be made. 

It appears worthwhile to study this problem with new efforts, because its solution will 
be toa large extent equivalent to the solution of Prar dtl’s boundary layer equation, but needs 
no restrictions as to the preseribed pressure distribution p(x), and will in principle be applicable 
to the more important case of turbulent boundary layers too. Since we want to include the 
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compressible case with heat transfer at the wall, we must simultaneously consider the temper- 
ature boundary layer. 

The infinite system of equations derived from Prandtl’s equation contains no terms which 
take into account an energy flux through the wall. Therefore, only the momentum equation 
v=0 and the mass-flow equation v=o of this system can be used in connection with heat 
transfer problems. All other equations are worthless in this respect. The essential step to 
get out of this difficulty is to include the partial differential equation for the temperature field 
in this problem by transforming it into an infinite system of ordinary differential equations 
similar to that obtained from Prandtl’s equation. 


3. Derivation of Integral Conditions From the Partial Differential Equation for 
the Temperature Field 


The general partial differential equation for the temperature field in a two-dimensional 
steady state flow may be written as follows, if the enthalpy ¢=e,T7 is introduced: 


ray) ay) dp ou\- O h Oo 
pu t+ pr su tu(—) + ~ | ‘ 
Ou Ov da Ov Oy \c, OY 
Here \ denotes the heat conductivity (keal/m sec " K), ¢, the specific heat at constant pressure, 

vu, ¢ must satisfy the continuity equation 


Olpu) , Apr) 
+ 


(). 
Ou Ov 


te 


Now, if we multiply eq (22) by 7, and eq (23) by ?*!/y+1 and add the multiplied equations, 
we obtain 


O(pui’*!) — Ol(pri?*!) dus . Ou \-. O dh O \. 
oo + — (y-+-1) —Upgs i? u( ) ie ( )i . 24 
ou Ov L dr Ov Ov \c, OY 
By partial integration of eq (24) with respect to y from y=0 to y=@ and by substitution for 


r from eq (23) we have 


d ( “+1 ee ; d | "yt ]) a) ‘ ; d ( ‘ i a . - 
— (pur gigi My i+ (psi i) 2 (pu ea) ae } ipslso} 
de a p p 5 Y dz p | y J, p f 7] 2 \f é 
1) ” du; | ~ nit . du ris | j r Oo; | , Ou ) h ae 
= (p-+ —pslls (i gis)? sls UstsO- ¢ Mi : a P -*) 
p dr, 0 - oe dr 0 f. Ov T oO: Y 


The terms with the factor 6 vanish because of the following known relation in gas dvnamies. 


| di; ue ] du; 


: +. 0), (26 
/ da is Us dar 


which follows after a short intermediate calculation from 


\J* - a*—critical sound velocity 28) 


1 dM; 1 du; ; 
Ms dx uy dr ( I+ Mj ): (29 


Using notations for the integral terms: 


J = ( , ; ) | ] u~ |e i ( ; ) | ] dy=% (30 


(@ can be replaced by 6, because the integrand is almost vanishing for y26) 


"'/ ps T u 
(2) Za {me 
Jo p Tua Us 


~ 
— 


‘ ( Ps ) O T Ol 7 T;) du 7 99 
| p ou r, Olu Ms Yes FS 4 
with 
Ou ou eC 
T= Oy" Tn | Al ’ 7 P. constant Prandtl number, (33 
and 
T / Ps Op 9 
a ae ( because Se 0 ): (34 
equation (25) may be rewritten as follows: 
‘ 
/o, dus | —M; k—1 T, y 
3 +9, : ~==—(y+1) —— Mj 2 Ze, + ] 7 . 39) 
dd dad 7 és su P . 
om P Mj 
L 2 


¢,, «and y, (with vy as an arbitrary number) are an infinite number of unknowns in this infinite 
system of ordinary differential equations. The real unknowns, however, are the distributions 
u(ry) u(r) and T(r,y)/Ts(2) which we are trying to determine by means of the infinite number 
of integral conditions as expressed in the system (35). Now, before going to the solution, let 


us transform it into a somewhat more convenient form. First a new variable is introduced: 


y= ?» ’ (36) 
; b.(k 1) Af? 
where 6» is the displacement thickness, defined bv 
"8 pu u . _ 
6 | —( ]— )dy to see (11)). (34) 
7 0 pss Us : 


It is to be remembered that », is finite when 1/4;->0, while ¢, has the disadvantage of becoming 


zero When \/;-+0. Furthermore, 


, 
- 


pd l 4 (k 1) ——, 6; | (1— eS Jey, (338) 
0» 0 


T ¥ 


Yr Ny ‘ 
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With those notations (35) can be rewritten as follows: 


V,=(y+1) 
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In principle a possible way to solve tne problem of determming the distribution u(r,y)/u;(x) 
and p(2,y)/ps(#) is to use polynomial expressions which fulfill the boundary conditions of the 
velocity and thermal boundary layer, by choosing 7 coefficients which satisfy 7 equations of 
the system (40). It is obvious that in the present case the procedure is more complicated 
than in the case of system (11). However, it may be possible to overcome by taking the solu- 
tion from the exact theory for the temperature profile 7/7; or the density profile p/p;= 73/7. 
Although the solution is exactly valid only for the special case of the flat plate (dp/dc=0), it 
may be considered as a good approximation even for the more general case with dp/de40 
(see for instance H. Schlichting [21, p. 282]). An essential property of this solution for 7/7; is 
that it depends not upon z but only upon w/v; with the Mach numbers 7; and the heat transfer 
coefficient ep as parameters. 

With reference to a paper of von Driest [23] this solution for 7/7; for the case dp/dr—0 


can be written as 


T : k—1 ° 
T. =1+-¢7(f,—1)+-7r > M?(1—f.) (41) 
with 
u u fu i ; 
h(P)=", f("P.)~(") for 0s<P,<1, (42) 
q Us Us ‘ U5 U5 
T.—T; 
r= pp —recovery factor, (43) 
T,—T, ‘ 
Cp=—,, “=heat transfer parameter. (44) 
7; 


Thus the temperature and density profiles occurring in the integral conditions (36) to (39) 
are uniquely correlated to the velocity profiles w/145, which now are the only unknown functions 
in the problem. So far all considerations and derivations are valid both for the laminar and 
turbulent case. For the solution of the laminar case, we have only to introduce an expression 
for the velocity profiles. The viscosity, as principal parometer of all boundary layer theories is 
known to be dependent only upon the temperature 7, for instance according to Sutherland's 


law. 


Lecce 
= a oF 


m 7 )" (49) 


(with C=113°K for air, 7, being a reference temperature) or to its more convenient approxi- 


mation 
m T\* “i 
Me ( T ) (46) 


’ 


l 
with 


w= — a] (47) 


| 
~ 
—~ 


(w~0.80 for temperatures 7° between 233° and 373°K, w=0.5 for very high temperatures, 
w=1.5 for very small temperatures near the zero point). 

In the case of turbulent boundary layers which are mainly considered in this report, an 
empirical relation for the effective viscosity (as the sum of the apparent and molecular viscosity) 


must be introduced besides the assumption for the velocity profiles. Before this is done, a 
simplification of the systems (11) and (40) will now be considered. 
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4. Special Case of a Two-Equation System 


4.1 General Considerations 


As a first step of evaluation of the possibilities involved in systems (11) and (40) we will 
make a so-called one-parametric representation for u/us, viz: 


= ( * H(z), (48) 


Us Us \6(4) 


where //(2) is the form parameter of the velocity profile (as one unknown of the problem) and 
6(r) is the boundary laver thickness (as second unkown). Therefore, to be correct, the relations 
(48) should be denoted as a two-parametric one. 

It must be pointed out that 6(7) may be replaced by another characteristic length of the 
boundary laver, preferably by the momentum-loss thickness 6,(7), thus avoiding any assumption 
as to the point where w/us reaches the value 1. 

Since we have now decided to use only two unknowns, //(7) and 6.(7), we need only two 
equations besides the empirical law for effective viscosity. 

We choose the simplest equations (namely those for y=0) from both the system (11) based 
upon Prandtl’s equation (1) and the system (40) based upon the general temperature field 
equation (2).° 

These equations are: 

a. The momentum integral condition from system (11). This is the only equation of the 
system (11) which, besides the mass-flow equation, is valid even in the case with heat transfer at 


the wall. 














- ‘ lo " uy 
Jor to 7 +25] - +é=—0 |> (49) 
0 Us 
with 
T - 
eé=— 9 (50) 
pauls 
. *8 ou u 
Jo | : | 1 iz b», (51) 
Jo Pals Us . 
°, pu . 
Mo | (1— )dy 5). (52) 
J pauls 


b. The energy integral condition from system (40) 


; us , Vo - 
Ny no? a - 0) (53) 
Us 0» 
with 
cs u p 
| ( 1— )dy 
] J Us Ps . (54) 
No = > - ’ ~ 
(k—1) M3 5, 
; , 6 : ms 
b=1+- (k4-1) Mj? —- (independent of »), (55) 
0» 
*) 
T, T u 1é a 
Som”. | d("")+m—3- }, (56) 
P55 J0 Tx Us -* 
5 The value of the simple mass-flow equation (19), which may also be obtained from system (40) by putting y= — ©, appears to be limited, but 


has not yet been studied in detail. 
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4 — 4) 
ah on rT; = VE 
s=P, l ( = ¥ a > )- is (for P,=0,8 to 1). (58) 


For the practical application of the two systems (49) and (53) it was found to be very 
helpful to evaluate once for all the universal functions 6,, 4), m, &, and Wo for the interesting 
range of the parameters /7, \/;, and @ (or er). Such an evaluation will be based upon a certain 
physically reasonable assumption about the velocity profile, as represented by (48). These 
assumptions for u/us; must of course be different in the laminar and turbulent cases. In the 


* us 


turbulent case, the empirical laws for the wall friction 7, and the dissipation integral rou 


must also be introduced. The evaluation for the laminar boundary layer contains no diffi- 
culties, and will be dealt with in a separate report. Only the more complicated case of the 
turbulent boundary laver, which needs some special considerations will be treated in this report. 


4.2. Two-Equation Method for Turbulent Boundary Layers 


We now consider: a) A suitable expression for empirical turbulent velocity profiles, 
b) an empirical law for the wall friction in a compressible turbulent boundary laver with heat 
transfer, ¢) an empirical law for the dissipation in a compressible turbulent boundary laver 
with heat transfer, and d) a method of solution of the svstem of equations (49) and (53 


a. The Turbulent Velocity Profile 


Detailed investigations of many authors (23, 24, 3] have shown that the turbulent boundary 
laver can be divided into several (usually into two) regions with respect to the direction normal 
to wall, in each one of which different laws for the influence of the outer pressure distribution 
and the Reynolds number are valid. Thus a general representation for the velocity profile 
requires a somewhat complicated formula. If only mean properties (integral quantities) of 
the boundary laver are considered and put into correlation as in our case, it will be sufficient to 
have an expression which contains less details but includes the whole range of the boundary. 

Such a formula, which is proved valid even at high Mach numbers, is 


where & is the form parameter with values approximately 0.1 and 0.7, and is a unique funetion 
of anv /7 to be chosen. 

A disadvantage of (59) is of course that at the wall y—0 the value (Ow Oy),~o is infinite for 
all 0k 1 and that furthermore, at y6—1, Ou Oy does not vanish as is physically required. 
But we have shown in previous publication |26, 27| that these facts do not essentially influence 
the integral terms of the approximate theory, since the wall-frietion law and the relation between 
the velocity profile u/u; and the shear stress profile 7/7, are separately introduced as semi- 
empirical laws. 

As an example of the suitability of this power law even for Mach numbers up to about 7 
and with heat transfer at the wall, we refer to figure 1 which is drawn from experiments of K. E. 
Wilson [25]. This simple power law is very adequate to the purpose of this approximate theory 
with regard to the possibility of closed evaluation of the integral terms. Figure 2 shows the 
order of magnitude of the errors in the ratio (0;/0)); plotted against (0; O.),, 
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Comparison of experiments (according to 7, Rotta and 


u 
hk. Weighardt) with approximation power law — 


(fs H)- 


caused by the approximation of the experimental velocity profiles with the power law calculated 
for the incompressible case. We see that the maximum error is about 1 percent. The neglect 
of a certain influence of the Reynolds number on the velocity profile causes errors less than 1 
percent. Therefore, the use of the power law (59) for our approximate theory appears to be 
justified, 

Some improvement at the outer edge region of the boundary laver can be reached without 
much complication by taking the following form: 


: =tanh |b log (14 ay)|, : (60) 
Us 
where a and 6 are two form-parameters, comparable to 7/7, and H1,. a may be used for the pur- 
poses of the adjustment of scales in the one-parametric theory, and the adjustment of wall 
friction in the two-parametric theory. The relation (60) avoids the use of 6 which is physically 
vague, and yields automatically an asymptotic behavior near the edge of the boundary layer. 
At the same time it prevents an infinite slope 0u/Oy at y=0. 

Investigations of the practical importance of expressions (60) (which of course, excludes 
the possibility of closed integrations) are in progress. 


b. The Wall Friction Law 


For this we may refer to the author’s publications dealing with the case of compressible 
turbulent boundary layers without heat transfer [26, 27|. We now have only to investigate to 
what extent the relations found there will be influenced by heat transfer at the wall. 

As a summary of these earlier investigations we can state that the laws for turbulent mixing 
motion found for the incompressible case, are valid in the compressible case, provided that the 
fluctuation velocity is less than the sound velocity, in other words, if the velocity fluctuations 

are not accompanied by density fluctuations. Now, since the turbulent fluctuation velocity 
in boundary lavers is found to be small, say 10 percent or 20 percent of the main velocity us 
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Us ° ~ , 
may be increased to about 10 or 5, respectively 
as Pe 


at the outer edge of the boundary laver, 1; 
until the fluetation motion becomes compressible in character. Therefore we conclude that up 
to about Mach number 1/= 10 or 5, respectively, the laws for incompressible fluctuation may 
be considered as valid to a first approximation. Only the flow medium properties u and p are 
changed. These, however, are taken into account by a suitably chosen Reynolds number. 

As to the effect of heat transfer at the wall, this will change only the flow medium properties 
and therefore can be included in the compressibility effect and in the Reynolds number, which 
itself plays a decisive role in the wall friction and dissipation law.’ 

The momentum law (v=0 from system (11)) contains the term = +, psu2? which is the 
ratio of the wall friction +7, to twice the stagnation pressure. This ratio can be rewritten 


as follows: 


ou ) | ° u/Us) O( UW / Us) ] 
T re Ov 7 Ov 6, u Oly (65) ;) u 65 (H V/ 4 ( [ “( | uf )ly (1 
— 9404 8.0). 0») . (1) 
u;J* 


Cr a pails. hs, (52) "Jo Me 


ao 


My, 


A definition of a Revnolds number 


psu 50> . 
R; (62) 


Ky 


is suggested which contains the viscosity uw, at the wall and, by means of the integral term 
6, an averaged inertia term pjujd.. Indeed the momentum integral condition (49) brings into 
relationship the wall friction and the averaged inertia forces. The pressure forces are not 
of interest in this connection. Any other definition of the Reynolds number appears to be not 
suitable in the present problem. In the incompressible case these detailed considerations are 
of course not necessary. ; 

In the numerator of eq (61) we have the dimensionless slope of u(y) at y—O0. As seale for 
the length y the value (62) jeompr. = (62)¢ as a wholly geometrical quantity, thus concentrating 
all influences of Mach number and heat transfer in the universal function 6,/(4),;. This fune- 
tion is theoretically given, when expressions for w/v; and p/p;, (48) and (41), are fixed. Next 


OU / Us) l ( ?. ) , 
= (3) 
O(y/ (82); Jul; \paus 


of eq (61) as the wall friction law for incompressible flow. For turbulent boundary lavers the 


the following term is considered: 


slope [O(u/u5)/04/62);J, depends itself on the Reynolds number 23, in such a manner that 

Olu U5) 

(y (65) u alll) T ° 

‘ . _ ( ;) (ti4) 

if 3 ts, pau; 
with »=0.268 (according to Ludwieg-Tillman [1S8]) and 
a( F1) =0.123-107 8° 1? (65) 

as an empirical function, which expresses the influence of shape of a 15. 


Therefore, we can finally write 


* A certain confirmation for these considerations may be seen in figure 1, which shows that the velocity distribution and therefore the related 


shear stress distribution r/r« too, is the same as in the incompressible case 


64 





tr, a@(H) 6 
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c. The Dissipation Law 


A generalization similar to that found for the wall friction can be derived for the dissipation 
in a compressible turbulent flow with heat transfer. The energy integral condition (53) con- 
tains the term 


| tdu=D (dissipation). (67) 


Written in dimensionless form, 


Dp f. "lg u Tw "log u 6» 
= d( ) ( = ) d( ) —- (6S) 
Ps pes Jo Te = \ Us pss 0 Te = \Us/_J (62); 


ta 


The term (7,,/psus"); is given with eq (64). The remaining integral term is found empirically 
by T. Rotta [33] and FE. Truckenbrodt |4] to be 
i u B( FT) 
d ( )= Py. (69) 
|, Te Us a(H) “** ' 


so that it finally vields 
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with 


D B(H) 0.0052 
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3 > 0. 168 0). 168 ( 71 ) 
sus), RY. I. 


as the part which represents the dissipation law for the incompressible flow. To be quite 
correct, a detailed consideration must be made on the definition of the Reynolds number 
used in eq (69). By definition the dissipation integral expresses the amount of energy dis- 
sipated by the viscosity into heat (with the intermediary of turbulence). The related Reynolds 
number should be more exactly defined with a mean viscosity 7 instead of uw», because the 
energy transfer by viscosity occurs within the whole boundary laver, although mostly near 
the wall. 
Therefore eq (69) should corres tly be written as 


ude IT) Pe a aT 
| T (7) ACD, Is. ( 


a | 
to 


‘wu 


when the definition (62) for the Reynolds number is used. This correction, however, was 
proved to be in general so small that it is not essential in the scope of this approximate theory. 


d. Simplified Representation for the Wall Friction and Dissipation Law 


The expressions for the incompressible wall friction and for the dissipation law were 
given by E. Truckenbrodt [4], who evaluated some semiempirical relations found by T. Rotta [3]. 
The results of the evaluation of E. Truckenbrodt are given in his two figures 3 and 4 
which we have reproduced in the present figures 3 and 4. These have been presented in 
order to investigate whether or not it will be possible to find a new suitable fitting for these 


65 














t 25us 
0 
“LT 7] | 
~ 
2.0 2 





1.0 ay 
ww 4 


Q 4 
ed . 














o.6}—} _— 
oa | “SE 
ES SS 2.0 














—rorTa '~ > 




















PY iy Bey fas | 0.2— D__o0076) 
2h | —"s ~~ ~ Aa RR? 
-0.665 ; 2 
Ze 00896-1072 9655/02) 
A uz Ree | R ost ; | + | dz, 
one oO 1 “de | oz, ie sor ? iF 
“107 2 5 /o* 2 5 105 
. . Fietre 4 Dissipation law for in- 
I IGURE 3 Wall triction law for in- compress hl, Kessiiaaitnael ai 
com presstbhle lurbulent houndar / — 
layers 
hick |] f J. Rotta-1 k 
Thiek line evaluation of e\periments hy 
I. Rotta-Truckenbrod Dashed line ippr mation proposed t 
Dashed lines evaluation of experiments o Trackentroed 
Luiwiee-Tillmann Dashed-dotted line: new approximati 
Dashed-dotted line: new approximation wit! 0.2 nower 1 
0.2 power law 


two empirical laws which permits the use (for the sake of simplification) of only one power 


of the Revnolds number instead of two (0.268 and 0.168). 
We find that if we take the power 0.2 in both cases the fitting of the experimental result 
in the most important range between /j.—10° and 10° will be rather satisfactory. The laws 


(66) and (70) may then be written as 


alae +) ) - 


rT, 0.0636 +10 “3 6b. indi 
pau; ie; * (6, ” 
DD O.OO7T61 6, =4 

pss ae )6| 6G 


This simplification prevents the universal function Wy of equation (56), from being dependent 
upon Rs, besides 77, MM; and 6. Thus an important simplification in the representation of 
this function seems possible without noticable loss of accuracy. Detailed investigations on 
the reliability of this simplification are now in progress. 

This simplification will be important if the use of higher) momentum equations 
(v=2,3,4, ... )is made. In this case all empirical parts of the wall friction and dissipation 


law are concentrated in the term r,,p v5, and the integrals of type 

adh 5 , T Ti . fa) r O 7 7 

| (f ) d( ) and | (* )v du 
0 p T. u Jo p Oul 7, OUuU 


are dependent only upon the shape parameter /7, the Mach number 1/5, and the heat transfer 
parameter #@ Therefore they can be evaluated, when the family of velocity profiles is chosen. 


66 


4.3. Evaluation of the Universal Functions of the Two-Equation Method for 
Laminar and Turbulent Boundary Layers 


The evaluation of the universal functions occurring in the two-equation method as described 
in section 4 has already been made in the author’s publications |26, 27] for the case of an insu- 
lating wall. 

The evaluation of these functions for the case with heat transfer, as suggested by the 
author, is now under way at the DVL Institute for Applied Mathematics and Mechanies, 
Freiburg Br, by considering first the turbulent boundary laver. 

A preliminary evaluation for this case with slide rule accuracy has been available for 
about one vear and was used for calculating some illustrating examples which will be discussed 
in the following section. 

This preliminary evaluation showed that the effect of heat transfer on all universal func- 
tions can be linearized in a rather good approximation, such that they may be represented 
in the form 


FV TM,,0) = FV T,M,0) +0-GUT, M,3). 75 


(4) 


This result will be checked after the completion of the final evaluation. 
4.4. Examples Flat Plate 


For the sake of completeness we present here again an important result obtained for 
the laminar case with insulated wall, figures 5 and 6. Pohlhausen’s method based only upon 
the momentum integral condition (49) and the wall condition (13) gives /7=1.572=const., 


for the ease of dp/dr—0 at each Mach number. The consideration of the energy integral 
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3) instead of eq (13) gives a variation of /7 with 17* and a corresponding change 
in the behavior of the friction coefficients ¢,/c,; (local at same values R;,) and Cp/Cp; (averaged) 
with the Mach number \/*. It is very interesting that the new behavior of e,/e, (MF) is 
almost coincident with the results of an exact theory of H. Wendt [8]. This result) shows 
the effectiveness of the energy integral condition. 


condition (: 


The related results for the compressible 
turbulent boundary layer with and without heat transfer are plotted in figures 7 and 8. 
The agreement with available experimental data is seen to be good. 
Figures 9 and 10 are concerned with the question of how the turbulent separation point 
ina retarded flow of the type 
un, M? 
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is influenced by the Mach number and by heat transfer. From figure 9 it is seen that the 
separation point (7//),.) is almost independent of variation in the Mach number between 
M*—0 to 2.40 (\7.=0 to 12). Figure 10 shows that by cooling the wall in such a way 
that the wall temperature 7, is kept equal to rs (- =1 ), the separation is completely prevented 
even for the high value MW" =2.40 (M.=12). On the contrary, heating of the wall moves 
the separation point upsiream. For the theroetical results of the figures 9 and 10, no experi- 
mental confirmation is vet available. 


5. Conclusions 


A survey of approximation theories for compressible laminar and turbulent boundary 
lavers with heat transfer at the wall and arbitrary pressure gradient in the flow direction 
shows that only a two-equation method has been used up to the present, although an infinite 
number of equations is available. A new derivation for the two infinite systems of integral 
conditions based upon the partial differential equations for the velocity and temperature 
field is given. 

An investigation of the significance of the equations with higher order shows that their 
property of deseribing typical boundary effects vanishes more and more with increasing number 
vy. Henee, it is not surprising that the results obtained from the two-equation method with 
v=0 in both the svstems (11) and (40) are good. However, the consideration of one or two 
equations more would be desirable, at least for a check of the results of the two-equation 
method in important special cases. 

Preliminary investigations have to show how the higher »-order integral terms containing 
the shear stress distribution 7/7, can be evaluated in the case of compressible turbulent 
boundary lavers with heat transfer. It has been found that this is possible if the power of 
the Reynolds number is chosen to be equal for both the wall friction and dissipation law. 
We have chosen 0.2. In spite of this simplification in these two empirical laws, they still fit 
the data rather well. Further work on this matter is in progress. 

This two-equation method now available with preliminary evaluated universal integral 
terms has been applied to the study of the dependence of the friction coefficient on compres- 
sibility and heat transfer. Good agreement with experimental data has been found. 
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A set of thirteen test matrices typical of those occurring in 
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inverses and various error estimates were computed. 
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On the convergence of the Rayleigh quotient itera- 
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putational advantages for more than three. 


Mode expansion in the low-frequency range for 
propagation through a curved stratified atmosphere, 
H. Bremmer, J. Research NBS 63D, No. 1,74 (1959). 
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ionospheric propagation at low frequencies. The complex 
problem dealing with two media, viz. a homogeneous earth 
and a surrounding stratified atmosphere, leads to intractable 
expressions. However, as the influence of the earth may 
be accounted for by an approximate boundary condition at 
the earth's surface, the problem is then reduced to that 
of the outer medium only. The coefficients of the mode 
expansion for this simplified problem will be derived while 
taking into account the earth’s curvature; however, the 
latter proves to be negligible under very general conditions, 
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